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^include <cstring.h> 
#include <:f stream. h> 
#include <:stdio.h> 
# include <:steiib.h> 
# include 'cmath.h> 

#de£ine TAB '\t' 
#define ARRAYHASH 256 



char base_f roTn_num (long num} 
{ 

switch (nuin%4) 
{ 



case 0 
case 1 
case 2 
case 3 



return ( 'A' ) 
return ( ' C ' ) 
return { * G ' ) 
return CT' } 



default: return ('A') 



long transform (string StTest} 
{ 

Test . to_upper ( ) ; 

for (long i=0; i<Test .length {) ; i++) 

{ 

if (Test [i]==*A' ) 
Test [i] = 'T' ; 

else 

if (Test [i] =='C' } 
Test [i] = 'G' ; 

else 

if {Test [iJ=='G' ) 
Test [i]='C' ; 

else 

if (Test [i] =='T' ) 

Test [i]='A' 

} 

return (TRUE) ; 



class ProbeSynth( 
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public: 

string Probe; 

string Original; 

string Name; 

atring Rename; 

long Location; 

long Unit; 

long Atom; 

// done with fancy 

char * Synth; 

long SynthLength; 

// shift by 4,8,12 

long SynthModif ier; 
long MutVal; // mutation position 

ProbeSynthO ; 
-ProbeSynth ( ) ; 
Destroy () ; 

Allocate (long) ; 
Synthesize {) ; 
SynthesizeFluff (long) ; 
SynthesizeFastFluf f (long) ; 
SynthesizeMutant (long) ; 
SynthesizeMisMatch(long} ; 
CompSynth ( ) ; 
char GetSynth(long) ; 

Set Synth {long, char) ; 
Zerocost ( ) ; 

Distance (ProbeSynth *) ; 
Output (of stream &) ; 
List {of stream &) ; 

}; 

ProbeSynth: : SynthesizeMutant (long MutateValue) 
{ 

long synthflag, done, fwdflag; 
long probecount, cyclecount, qloop; 
long start, finish; 
long re Synth; 

start = 0; 

finish = 0; 
Allocate (4 *Probe. length 0 +4) ; 

probecount = 0 ; 

this->MutVal = MutateValue; 

cyclecount = start; 

done - FALSE; 

synthflag = FALSE; 

fwdflag = TRUE; 

while (Idone) 
{ 

if { probecount ==Mutate Value) 
{ 

if (probecount <Probe. length (} -1) 

{ 

resynth = cyclecount; 

cyclecount=resynth+5; /*leave room for mutant*/ 
/♦synthesize mutant*/ 
synthflag = FALSE; 

for {qloop=re3ynth; qloop<cyclecount-l ; qloop++) 
{ 

if (base_from_num{ qloop) == Probe [probecount] ) 
{ 

// let us know that this is a mutant (and hence different) 
// >if< we were doing ACGT, filling in these with all 4 
// would let us match probes well, since geographic position 
// would matter. But we're only doing single mismatches, 
// for GE, so that doesn't work. 

SetSynth (qloop , ' * ' ) ; 

synthflag = TRUE; 

) 



cyclecount-- ; 

} 

else 
{ 



else 

Set Synth (qloop, '*'); 

} 

if (isynthflag) 

done = TRUE; 
probecount++; 



/*cancel out cycleco'unt++ later*/ 



if (base_from_num (cyclecount) == Probe [probeco\int3 ) 



Set Synth (cyclecount , base_from_nuin (cyclecount) ) , 
synthflag = TRUE; 
probecount ++ ; 

} 

} 

if (probecount >Probe . length (} -1) 
done = TRUE; 

if (finish>0 && synthflag) 
{ 

fwdflag - FALSE; 

probecount = Probe . length () - 1 ; 

cyclecount - finish; 

} 

else 
{ 

cyclecount ++; 
} 

1 

else 

if (base_from_num( cyclecount) ==Probe [probecount] ) 
{ 

if (fwdflag) 

{ 

SetSynth( cyclecount , base_from_nura( cyclecount) ) , 

probecount++; 

cyclecount++; 

if (probecount >= Probe - length ( ) ) 
done = TRUE; 

} 

else 



cyclecount-- ; 
} 

else 
{ 

cyclecount++; 
cyclecount-- ; 



Set Synth (cyclecount , base_f rom_nuin (cyclecount) ) 
probecount- - ; 

if {probecount <=MutateValue) 
done = TRUE; 



if (fwdflag) 
else 



this->SynthLength = cyclecount; // shorten 'search' space 
return (TRUE) ; 



ProbeSynth: : SyntheaizeMisMatch (long MutateValue) 



long synthflag, done, fwdflag; 
long probecount, cyclecount, qloop; 
long start, finish; 
long resynth; 
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start = 0; 
finish = 0; 

Allocate { 4 *Probe. length (} +4) ; 

probecount = 0 ; 

cyclecount = start; 
this->MutVal = MutateValue; 

done = FALSE; 

synthflag = FALSE; 

fwdflag = TRUE; 

while (idone) 

{ 



different) 



if ( probecount ==Mut ate Value) 
{ 

if (probecount<Probe. length (} -1) 
{ 

resynth = cyclecount ; 

cyclecount =resyiith+5; /*leave room for mutant*/ 
/*synthesi2e mutant*/ 
synthflag = FALSE; 

for (qloop=resynth; qloop<cyclecount-l ; qloop++) 

if {base_f rom_nuin{qloop) ==Probe [probecount] ) 
{ 

// let us know that this is a mutant (and hence 

Set Synth (qloop, '#'); 
synthflag = TRUE; 



cyclecount-- ; 



else 
{ 



if (! synthflag) 

done = TRUE; 
probecount ++; 

/♦cancel out cyclecount++ later*/ 



if (base_f rom_num{cyclecount) == Probe [probecount ] ) 
{ 

SetSynth {cyclecount, base_froni_nuni (cyclecount) ) ; 
synthflag = TRUE; 
probecount ++; 

} 

) 

if (probecount >Probe . length 0 -1) 
done = TRUE; 

if (finish>0 StSt synthflag) 
( 

fwdflag = FALSE; 

probecount = Probe . length 0 - 1 ; 

cyclecount = finish; 

} 

else 
{ 

cyclecount ++; 



else 



if (base_froin_num (cyclecount ) ==Probe [probecount] ) 
{ 

if (fwdflag) 
{ 

SetSynth (cyclecount, base_from_num (cyclecount) ) ; 
probecount ++; 

cyclecount ++; 

if (probecount>=Probe.length() ) 
done = TRUE; 

} 

else 
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cyclecount- - ; 



} 

else 



cyclecount++; 
} 



SetSynth (cyclecount , base_from_nuni (cyclecount) ) 
probecount-- ; 

if {probecount<=Mut ate Value) 
done = TRUE; 



if (fwdflag) 
else 



cyclecount-- ; 



} 

this->SynthLength = cyclecount; // shorten 'search' space 
return (TRUE) ; 

} 

char 

ProbeSynth: :GetSynth (long Which) 
{ 

if (Which < SynthLength && Which>-1) 
return (Synth [Which] ) ; 

; el se 

return { ' . ' ) ; 

} 

ij J 

== ProbeSynth: : SetSynth( long Which, char What) 

:z { 

if (Which>-1 && Which < SynthLength) 
;g Synth [Which] = What; 

return (TRUE) ; 

;i; > 

ProbeSynth: :ZeroCost ( ) 

;U { 

■■sJ return (Probe . length ( ) ) ; 

m } 

•~ ProbeSynth: : Distance {ProbeSynth *Destination) 

m { 

l"^: static char testchar; 

~ static long count; 

^5=1 static long minlen; 

static long i; 

count = 0 ; 

minlen = min (Destination- >SynthLength, this- >SynthLength) ; 

for (i=0; i<minlen; i++) 
{ 

testchar = De stinat ion -> Synth [i] ; // can guaremtee i>0 
if (testchar !='.') 

if (this->SynthIi] i=testchar) 

{ 

if (testchar == '#' ) 

count+=2; // heavily penalize mutants 

else 

count ++ ; 

} 

} 

for (; i<Destination->SynthLength; i++) 
{ 

testchar = De3tination->Synth [i] ; // can guarantee i>0 

if (testchar != ' . • ) 

{ 

if {testchar == '#') 

count+=2; // heavily penalize mutants 

else 
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count ++; 

} 



return (count) ; 



ProbeSynth: Synthesize 0 

{ 

long Size = CompSynthO ; 
Allocate {Size) ; 

long length, probeloop, synthloop; 
char testchar; 

length = Probe . length () ; 

probeloop = 0; 

synthloop = 0; 

while (probeloop-clength) 

{ 

testchar = base_froin_^nurrt (synthloop) ,- 

if ( test char==Probe [probeloop] ) 

{ 

Synth [synthloop] = Probe [probeloop] 
probeloop++; 

} 

else 

Synth [synthloop] = ' . ' ; 
synthloop++; 

} 

return (TRUE) ; 



ProbeSynth: : SynthesizeFluf f (long MutateValue) 
{ 

long length, probeloop, synthloop, testloop; 
char testchar; 

length = Probe . length ( ) ; 

long Size ~ 4*length+4; 

Allocate (Size) ; 
this->MutVal = MutateValue; 

probeloop = 0 ; 

synthloop = 0; 
for (probeloop=0 ; probeloop< length; probeloop++) 

{ 

synthloop= probeloop*4; 
testchar = Probe [probeloop] ; 
for (testloop = 0; test loop<4 ; testloop++, synthloop++) 
{ 

if (base_from_nuTn( synthloop) ==testchar) 



{ 
} 

} 



if (probeloop ! =MutateValue) 

Synth [synthloop] = testchar; 

else 

Synth [synthloop] = '#'; 
else 

Synth [synthloop] = ' . ' ; 



} 

return (TRUE) ; 

} 

ProbeSynth: : SynthesizeFastFluf f (long MutateValue) 

{ 

long length, probeloop, synthloop, testloop; 
char testchar; 

length = Probe . length () ; 
long Size = length+1; 
Allocate (Size) ; 
this->MutVal = MutateValue; 
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probeloop = 0; 
synthloop ~ 0; 

for (probeloop=0, synthloop = 0; probeloop< length; probeloop++, synthloop++) 
{ 

Synth [synthloop] = Probe [probeloop] ; 
if (probelocp==MutateValue) 
{ 

synthloop++; // account for mutant being twice as bad 
Synth [synthloop] =Probe [probeloop] ; 

1 

} 

return (TRUE) ; 



ProbeSynth: : Comp Synth () 
{ 

long length, probeloop, synthloop; 
char test char; 

length = Probe . length ( ) ; 

probeloop = 0 ; 

synthloop = 0; 

while (probeloop<length) 

{ 

testchar = base_frora_nuin (synthloop) ; 
if (testchar==Probe [probeloop] } 
{ 

probeloop++r 

} 

synthloop++; 

} 

return (synthloop) ; 



ProbeSynth .- ProbeSynth ( ) 

( 

■= Probe = 

Name = "None"; 

Location = 0; 
\_ Synth = NULL; 

ij3 SynthLength = 0; 

SynthModif ier = 0; 

'A ) 

ProbeSynth: : Allocate (long Size) 
{ 

Destroy () ; 
Synth = new char [Size] ; 
for (long i=0; i<Size; i++} 
Synth [i] = ' . ' ; 
SynthLength = Size; 

} 

ProbeSynth: :Destroy() 
{ 

if (Synth !=NULL) 



{ 



delete n Synth; 
Synth = NULL; 
SynthLength = 0 ; 



} 



ProbeSynth : : - ProbeSynth ( ) 
{ 

Destroy () ; 

} 

ProbeSynth: : Output (of stream &OutStream) 
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,sst 



OutStreatn <<: this- >Locat ion « TAB; 
Out St ream <<: this->Original << TAB; 
OutStream <<; this->Nanie << TAB; 
OutStream << this-sRename « TAB; 
OutStreain « this->Unit << TAB; 
OutStream << this->Atom << endl; 



ProbeSynth: :List (of stream tOutStream) 
{ 

OutStream c< thi9->0riginal « endl; 

} 

class ProbeNode{ 
public : 

ProbeSynth *DataPointer ; 
ProbeNode * Previous; 
ProbeNode *Next; 

ProbeNode **PClosest; 
long plength; 
long marker; 
long Next Co St ; 
PointToData( ProbeSynth *) ; 
ProbeNode ( ) ; 
-ProbeNode ( ) ; 
Destroy 0 ; 
DestroyPClosest () ; 

DestroyData () ; 
AllocatePClosest (long) ; 

LinkPrevious {ProbeNode *); 
' ^= LinkNext ( ProbeNode * ) ; 

;0 I nsertNext (ProbeNode *) ; 

Initialize (} ; 
'.l". ZeroCost () ; 

I y Copy ( ProbeNode * ) ; 

Distance {ProbeNode *); 

i3 

\J} ProbeNode: : Copy (ProbeNode *Original) 

;„ if (Original ! =NULL) 

DataPointer = Original ->DataPointer; 

i3 else 

return (FALSE) ; 

// note - do not copy Pclosest, plength - don't apply to copies 
// do not copy previous & next, because they won't correspond. 



ProbeNode: : PointToData (ProbeSynth *Data} 
{ 

DataPointer = Data; 



ProbeNode : : ZeroCost ( ) 
{ 

return (DataPointer- >ZeroCost () ) ; 

1 

ProbeNode: : Distance (ProbeNode *Destination) 
{ 

return (DataPointer- >Distance (Destination- >DataPointer) ) , 



ProbeNode : : ProbeNode ( ) 
{ 

DataPointer - NULL; 
Previous = NULL; 
Next = NULL; 

PClosest = NULL; 
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marker = 0; 

plength =0; 
NextCost = 0; 

} 

ProbeNode : .-Destroy ( ) 
{ 

DataPointer = NULL; 
Previous = NULL; 
Next - NULL; 

) 

ProbeNode: .DestroyPClosest () 
{ 

if (PClosest !=NULL) 

{ 

delete [] PClosest; 
PClosest = NULL; 
plength = 0; 

} 



ProbeNode : : DestroyData ( ) 
{ 

// dangerous - must be last called 
if (DataPointer !=NULL} 
{ 

delete [] DataPointer; 
DataPointer = NULL; 
} 

} 

ProbeNode: : AllocatePClosest (long Size) 
{ 

DestroyPClosest 0 ; 

PClosest = new ProbeNode * [Size] ; 
for (long i=0; i<Size; 

PClosest ti] = NULL; 
return (TRUE) ; 

} 

ProbeNode : : -ProbeNode ( } 
{ 

Destroy 0 ; 

} 

ProbeNode: : LinkPrevious (ProbeNode *Link) 
{ 

Previous = Link; 

} 

ProbeNode : : LinkNext ( ProbeNode *Link) 
{ 

Next = Link; 

} 

ProbeNode : : InsertNext {ProbeNode *NewNode) 
{ 

ProbeNode *Link; 

Link = Next ; 

Next = NewNode; 
NewNode-> Previous = Link->Previous; 
Next->Previou3 = NewNode; 
NewNode- >Next = Link; 



ProbeNode: : Initialize {) 
{ 

Previous = this; 
Next = this; 



} 



class TourClass{ 
public : 

ProbeNode *DataList; 
long Cost; 

TourClass () ; 

InitializeTour (ProbeNode *) ; 
DeleteCurrent ( ) ; 
UnlinkCurrent (} ; 
DestroyList () ; 

Insert Aftercurrent (ProbeNode *) ; 
long TestBasicInsertion (ProbeNode *) ; 
Rotate 0 ; 

Duplicate (TourClass &) ; 
InsertLeastCost {ProbeNode *} ,- 
InsertLeastCostFrotnPooI (TourClass &) ; 
InaertLeastCostByMarker (TourClass &) ; 
InsertLeastCostWithModifier (ProbeNode long) ; 
QuicklnsertLeastCost (ProbeNode *, long); 
-TourClass 0 ; 
Output (of stream &) ; 
List (of stream &) ; 

}; 

TourClass: : TourClass () 
{ 

DataList = NULL; 
Cost = 0; 

) 

TourClass: :Output (of stream &OutStream) 
{ 

ProbeNode *Start = this- >DataList ; 

this- >DataList- >DataPointer->Output {OutStream) ; 
Rotate 0 ; 

while (Start != this->DataList) 
{ 

this- >DataList ~ >DataPointer- >Output (OutStream) 
Rotate ( ) ; 

} 

} 

TourClass: ; List {of stream &OutStream} 
{ 

ProbeNode *Start = this- >DataList ; 

this->DataList->DataPointer->List (OutStream) ; 
Rotate { ) ; 

while (Start != this- >DataList ) 
{ 

this->DataList->DataPointer->List (OutStream) ; 
Rotate ( ) ; 

} 

) 

TourClass : : Duplicate (TourClass &TestTour) 
{ 

ProbeNode * Duplicate; 
ProbeNode *Start; 

this->DestroyList () ; 

Start = TestTour .DataList ; 
Duplicate = new ProbeNode; 
Duplicate->Copy (TestTour .DataList) ; 

InitializeTour (Duplicate) ; 
TestTour . Rotate ( } ; 

while (Start !=TestTour. DataList) 
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{ 

Duplicate = new ProbeNode,- 
Duplicate- >Copy (TestTour .DataList) ; 
InsertAfterCurrent (Duplicate) ; 
TestTour . Rotate { ) ,- 

Rotate < ) ; 

} 

return (TRUE) ; 

} 

TourClass: : insertAfterCurrent {ProbeNode *NewNode) 
{ 

ProbeNode *Link; 

if (DataList==NULL) 
{ 

InitializeTour (NewNode) ; 
return (TRUE) ; 

} 

Link = DataList- >Next ; 

DataList ->Next = NewNode; 

Link->Previous = NewNode; 

NewNode ->Next = Link; 

NewNode- > Previous = DataList; 

Cost = Cost - DataList- >NextCost ; 

DataList ->NextCost = DataList- >Distance (NewNode) ; 

Cost = Cost + DataList ->Next Cost ; 

NewNode- >Next Cost = NewNode- >Di stance (Link) ; 

Cost = Cost + NewNode ->Next Cost ; 

} 

long 

TourClass: :TestBasicInsert ion (ProbeNode *NewNode) 
{ 

long TestCost ; 

TestCost = Cost - DataList ->Next Cos t ; 

TestCost += DataList ->Distance (NewNode) ; 
TestCost += NewNode->Distance (DataList - >Next ) ; 
return (TestCost) ; 

} 

TourClass : : Rotate ( ) 
{ 

DataList = DataList - >Next ; 

} 

TourClass: ; insertLeastCost (ProbeNode *NewNode) 
{ 

ProbeNode *Start; 
ProbeNode *BestPlace; 
long TestCost, BestCost; 

Start = DataList; 

BestPlace = DataList; 
TestCost = TestBasicInsertion (NewNode) ; 

BestCost = TestCost; 

Rotate 0 ; 
while (DataList !=Start) 

{ 

TestCost = TestBasicInsertion (NewNode) ; 
if (TestCost<BestCost) 
{ 

BestCost = TestCost; 
BestPlace = DataList; 

} 

Rotate 0 ; 

} 

DataList = BestPlace; 
InsertAfterCurrent (NewNode) ; 

} 
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TourClass: : InsertLeastCostByMarker (TourClass &Source) 
{ 

// searches Source for undone ones, 
// adds the closest one to the new tour 

ProbeNode *Start; 
ProbeNode *BestPlace, *BestAdd; 
ProbeNode *NewNode, *Done; 

long TestCost, BestCost; 



Start = DataList; 

BestPlace = DataList; 
NewNode = Source .DataList ; 

Done = Source. DataList ; 
TestCost = TestBasicInsertion (NewKode) ; 
BestCost = TestCost; 
Rotate ( ) ; 

Source . Rot at e () ; 
long totalmarker = 0; 

while {Source . DataList ! =Done) ; 

{ 

NewNode = Source .DataList ; 
while (DataList ! =Start && ! NewNode -^marker) 
{ 

TestCost = TestBasicInsertion (NewNode) , 
=r" if (TestCost-cBestCost) 

BestAdd = NewNode; 
^ BestCost = TestCost; 

£ BestPlace = DataList; 

~ Rotate 0; 

y totalmarker = 1; 

a : Source . Rotate ( ) ; 

^ } 

; NewNode = new ProbeNode; 

*=i NewNode ->Copy (BestAdd) ; 

^ BestAdd- >marker = 1; 
^ DataList = BestPlace; 

InsertAf terCurrent (NewNode) ; 
~ return (totalmarker) ; 

i: } 

=t TourClass :: InsertLeastCostFroraPool (TourClass &Source) 
t 

// searches Source for undone ones, 
// adds the closest one to the new tour 

ProbeNode * St art; 
ProbeNode *Best Place, *BestAdd; 
ProbeNode *NewNode, *Done; 
long TestCost, BestCost; 



BestAdd = NULL; 

BestPlace = NULL; 

NewNode = Source . DataList ; 

Done = Source. DataList; 

TestCost = 1000; 

BestCost = TestCost; 

Source . Rotate ( ) ; 

long totalmarker = 0 ; 

while (Source. DataList ! =Done) 

{ 

NewKode = Source .DataList ; 

//cout << NewNode- >DataPointer->Probe << TAB << NewNode - >marker << endl 

if (! NewNode- >marker) 

{ 

TestCost = TestBasicInsertion (NewNode) ; 
if (TestCost <Best Cost J 

{ 

BestAdd = NewNode; 
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BestCost = TestCost; 

) 

totalmarker = 1; 

} 

Source . Rotate ( } ; 

} 

if (totalmarker && Best Add) 
{ 

NewNode = new ProbeNode; 

NewNode - >Copy ( Bes t Add > ; 

BestAdd->Tnarker = l; 

InsertAfterCurrent (NewNode) ; 

cout << NewNode- >DataPointer->Probe << endl; 

Rotate {); // go to NewNode as the favorite 

} 

return (totalmarker) ; 



TourClass; :lnsertLeastCostWithModifier( ProbeNode *NewNode, long ModMax) 

{ 

ProbeNode *Start; 
ProbeNode *BestPlace; 
long BestModif ier; 
long TestCost, BestCost; 

Start = Dataliist; 
BestPlace = DataList; 

TestCost = TestBasicInsert ion (NewNode) ; 
BestCost = TestCost; 
BestModif ier = 0; 

for (long ModLoop = 0; ModLoop< ModMax; ModLoop++) 
{ 

NewNode- >DataPointer->SynthModi£ ier = ModLoop ; 
Start = DataList; 
Rotate 0 ; 

while (DataList !=Start) 
{ 

TestCost = TestBasicInsert ion (NewNode) ; 

if {TestCost<BestCost) 

{ 

BestCost = TestCost; 
BestPlace = DataList; 
BestModif ier = ModLoop; 

) 

Rotate () ; 

} 

} 

DataList = BestPlace; 

NewNode->DataPointer->SynthModif ier = BestModif ier; 
InsertAfterCurrent (NewNode) ; 



TourClass QuicklnsertLeastCost {ProbeNode *NewNode, long SearchLevel) 
{ 

ProbeNode * St art; 
ProbeNode *Best Place; 
long TestCost, BestCost; 

Start = DataList; 
BestPlace = DataList; 

TestCost = TestBasicInsert ion (NewNode) ; 
BestCost = TestCost; 

Start = DataList; 

Rotate 0 ; 

long counter = 0; 

NewNode~>DataPointer->SynchModif ier * 0; 
while (DataList !=Start counter<SearchLevel) 
{ 

TestCost = TestBasicInsert ion (NewNode) ; 
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if (TestCost<BestCost) 
{ 

BeatCost = TestCost; 
BeatPlace = DataList; 

) 

Rotate 0 ; 
count er++; 

) 

DataList = BestPlace; 
InserCAf terCurrent (NewNode) ; 

1 

TourClass: .-DeleteCurrent () 
{ 

ProbeNode *Link; 
ProbeNode *Del; 

Link = DataList ->Next; 
Del = DataList; 

Cost = Cost - DataList ->NextCost; 

Cost = Cost - DataList ->Previous->Next Cost; 
if (Link==DataList ) 
{ 

DataList = NULL; 
delete Del; 

Cost = 0; 
return (TRUE) ; 

} 

Link->Previous = DataList ->Previous; 
Link->Previous->Kext = DataList- >Next ; 
DataList = Link; 
delete Del; 

Link->Previous->NextCost = Link- >Previous->Distance (Link) ; 
Cost = Cost + Link->Previous->NextCost ; 
return (TRUE) ; 



TourClass : :UnlinkCurrent () 

{ 

ProbeNode *Link; 
ProbeNode *Del ; 

Link = DataList ->Next; 
Del = DataList; 

Cost = Cost - DataList- >NextCost ; 

Cost = Cost - DataList- >Previous->NextCost ; 
if (Link==DataList) 
( 

DataList = NULL; 

Cost = 0; 
return (TRUE) ; 

} 

Link->Previous = DataList -> Previous ; 

Link->Previous->Next = DataList- >Next ; 
DataList = Link; 

Link->Previous->NextCost = Link->Previous->Distance (Link) 
Cost = Cost + Link->PreviQus- >NextCost ; 

return (TRUE) ; 



TourClass: : initializeTour (ProbeNode *Test} 
{ 

DestroyList () ; 
DataList = Test; 
Test->Initialize 0 ; 
Cost = Test->ZeroCo3t 0 ; 
Test->NextCost = 0; 
return (TRUE) ; 



TourClass: -.DestroyList ( ) 
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while (DataList J=NULL) 
DeleteCurrent {) ; 
) 



TourClass: : -TourClass ( } 
{ 

DestroyList {} ; 

} 

class TSPClas3{ 
public: 

TourClass DataSet; 
TourClass BestTour; 
TourClass CurrentTour; 
TSPClass 0 ; 
-TSPClass 0 ; 
LoadData (string) ; 
LoadMutantDat a (string, long) ; 
LoadMisMatchData (string, long); 
LoadMisMatchFluff Data (string, long) ; 
LoadMisMatchFastFluff Data (string, long) ; 
LoadExpressionData (string, long); 
LoadSingleExpressionData (string) ; 
LoadExpressionDataByUnit (string, long, long) 
LoadExpressionFluff Data (string, long) ; 
LoadChipData( string) ; 
GenerateTourBylnsertion 0 ; 
GenerateTourBylnsertionAndDeletion ( ) ; 
GenerateTourByClosestlnsertion ( } ; 
GenerateTourByClosestPool () ; 
GenerateTourBylnsertionWithModif ier (long) ; 
GenerateQuickTourBylnsert ion (long) ; 
ImproveTourByReplacement (long) ; 
OutputTour (string) ; 
AppendTour (string) ; 
ListTour (string) ; 
DestroyData 0 ; 
// Geni(int); 

}; 

TSPClasS: :TSPClass() 

{ 

} 

TSPClass: : DestroyData 0 
{ 

ProbeNode *DataStart; 

DataStart = DataSet .DataList ; 
DataSet .DataList ->DestroyData 0 ; 
DataSet .Rotate () ; 

while (DataStart ! =DataSet . DataList) 
{ 

DataSet .DataList->DestroyData( ) ; 
DataSet .Rotate 0 

) 

return (TRUE) ; 

} 

TSPClass : : -TSPClass ( ) 
{ 

//DestroyData () ; 

} 

TSPClass: : GenerateTourBylnsertion {) 
{ 

ProbeNode *DataStart; 
ProbeNode *TenipData; 
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DataStart = DataSet .DataList ; 
TempData = new ProbeNode; 
TeTnpData->Copy (DataSet .DataList) ; 
BestTour .InitializeTour (TempData) ; 
DataSet - Rotate 0 

long counter = 0; 

while (DataStart ! =DataSet . DataList ) 
{ 

TempData - new ProbeNode; 
TempData- >Copy (DataSet .DataList ) ; 

//cout « DataSet .DataList ->DataPointer->Probe « TAB << counter << TAB; 
BestTour . InsertLeastCost (TempData) ; 
//cout << BestTour . Cost << endl; 
DataSet .Rotate {) ; 
if (counter%100==0) 

cout << counter « endl; 
counter++; 

} 

return (TRUE) ; 

} 



TSPClass : iGenerateTourBylnsertionAndDeletion 0 
{ 

ProbeNode *DataStart; 
ProbeNode *TempData; 



DataStart = DataSet .DataList ; 
if (DataSet .DataList ->marker>0) 

DataSet .Rotate {) ; 
while (DataSet .DataList- >marker>0 DataSet .DataList !=DataStart) 
{ 

DataSet .Rotate 0 ; 

} 

if (DataSet .DataList- >marker>0) 

return (FALSE) ; 
DataStart = DataSet .DataList ; 
DataStart- >marker = 1; 

TempData = new ProbeNode; 
TempData- =>Copy (DataSet - DataList) ; 
BestTour. InitializeTour (TempData) ; 

long Unit = TempData->DataPointer->Unit ; 

DataSet -Rotate () ; 

long counter = 0 ; 

while (DataStart ! =DataSet . DataList ) 
{ 

if (DataSet .DataList ->DataPointer->Unit==Unit DataSet .DataList ->marker<l) 
{ 

TempData = new ProbeNode; 
TempData- >Copy (DataSet . DataList) ; 

cout « DataSet .DataList- >DataPointer->Name << TAB << DataSet .DataList- 
>DataPointer->Unit << TAB « counter « endl; 

BestTour . InsertLeastCost (TempData) ; 
DataSet .DataList- >marker = 1; 
//cout << BestTour .Cost << endl; 
} 

DataSet . Rotate ( ) ; 
count er++ ; 

if (counter%1000==0) 

cout << coxinter << endl; 

} 

return (TRUE) ; 

} 
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TSPClass : ;GenerateTourByClosestln8ertion() 
{ 

ProbeNode *DataStart; 
ProbeNode *TempData; 

DataStart = DataSet .DataList ; 
TetnpData = new ProbeNode ; 

TempData->Copy C DataSet .DataList) ; 
DataSet .DataList ->Tnarker = 1; // set on this course 

BestTour . iniCializeTour (TempData) ; 

long notdone = 1; 
long counter = 0; 
while (notdone) 
{ 

notdone = BestTour .insertLeastCostByMarker (DataSet) ; 
if (counter%100==0) 

cout << counter << endl; 
counter++; 

} 

return (TRUE) ; 

} 

TSPClass ; :GenerateTourByC lose St Pool () 

ProbeNode *DataStart; 
ProbeNode *TeTnpData; 

DataStart ~ DataSet . DataList ; 
TempData = new ProbeNode; 
TetnpData- >Copy (DataSet .DataList) ; 

DataSet .DataList ->marker = 1; // set on this course 
BestTour. InitializeTour (TempData) ; 

long notdone = 1; 
long counter = 0; 
while (notdone) 
{ 

notdone = BestTour . InsertLeastCostFromPool (DataSet ) ; 
counter ++; 

} 

return (TRUE) 



TSPClass: :GenerateTourSyInsertionWithModif ier (long ModMax) 
{ 

ProbeNode *DataStart; 
ProbeNode * TempData ; 

DataStart = DataSet . DataList ; 
TempData = new ProbeNode ; 
TempData- >Copy (DataSet .DataList) ; 
BestTour. InitializeTour (TempData) ; 
DataSet . Rotate ( ) ; 

long counter = 0; 

while (DataStart i =DataSet .DataList) 
{ 

TempData = new ProbeNode; 
TempData- >Copy (DataSet .DataList) ; 

BestTour . InsertLeastCostWithModif ier (TempData, ModMax) ; 

if {counter%100==0) 

{ 

cout << DataSet .DataList ->DataPointer->Probe << tab << counter << TAB; 
cout « BestTour .Cost << endl; 

} 

DataSet .Rotate {) ; 
counter++; 

} 

return (TRUE) ; 

} 
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TSPClass : : GenerateQuickTourBylnsert ion ( long SearchLevel ) 
{ 

ProbeNode *DataStart; 
probeNode *TeTnpData; 

DataStart = DataSet .DataList ; 
TempData = new ProbeNode; 
TempData->Copy (DataSet .DataList) ; 
BestTour . Initial izeTour (TempData) ; 
DataSet. Rotate 0 ; 

long counter = 0; 

while {DataStart !=DataSet .DataList) 
{ 

TempData = new ProbeNode; 
TempData- >Copy (DataSet .DataList) ; 

BestTour .QuicklnsertLeastCost (TempData, SearchLevel) ; 
if (counter%100==0) 
{ 

cout << DataSet .DataList ->DataPoinCer->Probe « TAB << counter << TAB; 
cout << BestTour .Cost << endl; 

) 

DataSet .Rotate () ; 
counter++; 

} 

return (TRUE) ; 



TSPClass : : ItnproveTourByReplacenient (long ReplaceSize) 
{ 

ProbeNode *DataStart; 
ProbeNode *TempData; 



long counter = 0 ; 

while (counter<ReplaceSize) 
{ 

TempData = BestTour . DataList ; 
if (TempData- >mar)cer<l) 

{ 

TempData- >marker = 1; 
BestTour .UnlinkCur rent () ; 
cout <<: TempData -'>DataPointer-=.Probe << TAB <:< counter << TAB; 
BestTour .insertLeastCost (TempData) ; 
cout << BestTour .Cost << endl; 
} 

else 

BestTour . Rotate ( ) ,- 

counter++; 
) 

return (TRUE) ; 



TSPClass: :LoadData( string FileName) 
{ 

// build the basic datalist 
if stream ExpStream; 

ExpSt ream. open { FileName. c_str () , ios; ;in) ; 

if { ExpSt ream. bad 0 ) 
{ 

cout « endl « "Experimental data file specified does not exist!"; 
cout << endl << endl; 
exit (0) ; 

} 

long probeloop = 0; 
string TestString; 
ProbeNode *NewNode; 

float test ; / / soak up number 

string junkgtring; 



{1 




18 



while ( lExpStream.eof {} ) 
{ 

ExpStreatn >> TestString >> junkstring » test; 
if (TestString. length {)>1) 

{ 

NewNode = new ProbeNode; 
NewNode->DataPointer = new ProbeSynth; 
NewNode- >DataPointer->Original * TestString; 
NewNode- >DataPointer->Probe = TestString; 
NewNode- >DataPointer- ^Synthesize {} ; 
this->DataSet . InsertAf terCurrent (NewNode) ; 
//cout << TestString c< endl; 

} 

} 

ExpStream . close ( ) ; 

cout << DataSet.Cost << endl; 

} 

TSPClass: :LoadMutantData( string FileName, long MutateValue) 

{ 

// build the basic dataiist 
if St ream ExpStream; 

ExpStream. open ( FileName. c_str() , ios : :in) ; 

if { ExpStr earn. bad (} ) 
{ 

cout <«: endl << "Experimental data file specified does not exist 
cout <<: endl << endl; 
exit (0) ; 

} 

long probe loop = 0; 

string TestString; 

ProbeNode *NewNode; 

float test; // soak up number 

while ( ! ExpStream . eof ( ) ) 
{ 

ExpStream >> TestString >> test; 
if (TestString. length 0 >1} 

{ 

NewNode = new ProbeNode; 
NewNode- >DataPointer = new ProbeSynth; 

NewNode- >DataPointer- >Probe = TestString; 
NewNode- >DataPointer->SynthesizeMutant (MutateValue) ; 
this->DataSet . InsertAf terCurrent (NewNode) ; 
} 

} 

ExpStream. close (} ; 
cout << DataSet.Cost << endl; 

} 

TSPClass: :LoadMisMatchData (string FileName, long MutateValue) 
{ 

// build the basic dataiist 
ifstream ExpStream; 

ExpStream. open (FileName. c_str{) , ios: tin) ; 

if (ExpStream. bad 0 ) 
{ 

cout << endl << "Experimental data file specified does not exist 
cout « endl << endl; 
exit (0) ; 

} 

long probeloop = 0; 
string TestString; 

ProbeNode * NewNode; 

float test; // soak up number 

while (! ExpStream. eof () ) 
{ 
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ExpStream >> TestString » test; 

if (TestString. length 0 >1) 

{ 

NewNode = new ProbeNode; 
NewNode->DataPointer = new ProbeSynth; 

NewNode- >DataPointer->Probe ~ TestString; 

NewNode- >DataPointer- >SyTithesizeMiaMatch (MutateValue) ; 

this- >DataSet -Insert Aftercurrent (NewNode) ; 

} 

} 

ExpStream . close ( ) ; 
cout << DataSet.Cost « endl; 

} 

TSPClass: :LoadMisMatchFluff Data (string FileName, long MutateValue) 
{ 

// build the basic datalist 
if stream ExpStream; 

ExpStream. open (FileName. c_str() , ios: : in) ; 

if (ExpStream.badO ) 
{ 

cout << endl « "Experimental data file specified does not exist 
cout << endl << endl; 
exit (0) ; 

} 

long probe loop = 0; 
string TestString; 

ProbeNode *NewNode; 

float test; // soak up number 

while (! ExpStream. eof {) ) 
{ 

ExpStream >> TestString >> test; 
if (TestString . length ( ) >1) 
{ 

NewNode = new ProbeNode; 

NewNode- >DataPointer = new ProbeSynth; 
NewNode- >DataPointer-:>Probe = TestString; 
NewNode->DataPointer->SynthesizeFluf f (MutateValue) ; 
this- >DataSet . InsertAf terCurrent (NewNode) ; 

} 

} 

ExpStream. close 0 ; 
cout << DataSet.Cost << endl; 

} 

TSPClass: :LoadMisMatchFastFluff Data (string FileName, long MutateValue) 
{ 

// build the basic datalist 
if stream ExpStream; 

ExpStream. open (FileName. c_str(} , ios: :in) ; 

if (ExpStream.bad( ) ) 
{ 

cout << endl << "Experimental data file specified does not exist 
cout << endl «:< endl; 
exit (0) ; 

) 

long probeloop = 0; 
string TestString; 

ProbeNode *NewNode; 

float teat; // soak up number 

while (! ExpStream. eof () ) 

{ 

ExpStream >> TestString >> test; 
if (TestString. length 0 >1) 
{ 

NewNode = new ProbeNode; 

NewNode- :>DataPointer = new ProbeSynth; 
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NewNode->DataPointer->Probe = TestString; 
NewNode->DataPointer->SynthesizeFastFluf f (MutateValue) ; 
this->DataSet . inaertAf terCurrent (NewNode) ; 

) 

} 

ExpStream. close 0 ; 

cout << DataSet.Cost << endl; 



TSPClass :; LoadExpresaionData (string FileName, long MutateValue) 
{ 

// build the basic datalist 
ifstreara ExpStream; 

ExpStream. open (FileName. c_3tr () , ios: :in) ; 

if (ExpStream.badO ) 
{ 

cout << endl << "Experimental data file specified does not exist!"; 
cout << endl « endl; 
exit (0} ; 

} 

long probeloop = 0; 

string TestString; 

ProbeNode *NewNode; 

float test; // soak up number 

long facevalue; 

long Unit, Atom; 

string Name, Rename; 

string Junk; 

long counter = 0; 

ExpStream >> Junk >> Junk >> Junk >> Junk >> Junk >> Junk; 

while (! ExpStream. eof 0 ) 

{ 

ExpStreani >> facevalue >> TestString » Name >> Rename » Unit >> At 

if (TestString. length C } >1) 

{ 

NewNode = new ProbeNode; 
NewNode- >DataPointer = new ProbeSynth; 
NewNode- >DataPointer->Original = TestString; 
transform (TestString) ; 

NewNode- >DataPointer->Probe = TestString; 
NewNode- >DataPointer->Naine = Name; 
NewNode- >DataPointer->Rename = Rename; 
NewNode- >DataPointer- >Location = facevalue; 
NewNode- >DataPointer->Unit = Unit; 
NewNode- >DataPointer->Atom = Atom; 

NewNode- >DataPointer- >SynthesizeMisMatch (MutateValue) ; 
this- >DataSet . InsertAf terCurrent (NewNode) ; 

} 

counter++; 

if (counter%100==0) 

cout << counter « TAB << Name << TAB << TestString << endl; 

) 

ExpStream. close ( ) ; 

cout << DataSet.Cost << endl; 

} 

TSPClass: ;LoadSingleExpressionData( string FileName) 
{ 

// build the basic datalist 
if stream ExpStream; 

ExpStream. open (FileName. c_str() , ios: --in) ; 

if {ExpStream.badO) 
{ 

cout << endl « "Experimental data file specified does not exist!"; 
cout « endl « endl; 
exit (0) ; 

} 

int probeloop = 0; 
string TestString; 
ProbeNode *KewNode; 
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float test; // soak up number 
int facevalue; 
int Unit, Atom; 
string Name; 
string Junk; 
int counter = 0; 

ExpStream >> Junk » Junk » Junk >> Jvmk » Junk » Junk; 
while ( lExpStream.eof (} } 
{ 

ExpStream >> facevalue » TestString >> Name >> Jimk >> Unit >> Atom; 
if (TestString. length {)>!) 
{ 

NewWode = new ProbeNode; 
NewNode->DataPointer ^ new probeSynth; 
NewNode->DataPointer->Original = TestString; 
transform (TestString) ; 

NewNode->DataPointer->Probe = TestString; 
NewNode->DataPointer->Name = Name; 
NewNode- >DataPointer->Location = facevalue; 
NewNode->DataPointer->Unit = Unit; 
NewNode- >DataPointer->Atom = Atom; 
NewNode ->DataPointer->Synthesize ( ) ; 
this->DataSet . InsertAf terCurrent (NewNode) ; 

} 

counter++; 
if (counter%100==0) 

n cout << counter << TAB << Name << TAB << TestString << endl ; 

ExpStream . close ( ) ; 
fi= cout << DataSet.Cost << endl; 

□ ) 

y TSPClass :: LoadExpressionDataByUnit (string FileName, long MutateValue, long UnitLimit) 
P { 

=5: // build the basic datalist 

~ ifstream ExpStream; 

ExpStream, open (FileName .c_str () , ios::in); 



if ( ExpStream, bad {) ) 
{ 

cout << endl << "Experimental data file specified does not exist!"; 
cout << endl <<; endl; 
exit (0) ; 

) 

long probe loop = 0; 

string TestString; 

ProbeNode *NewNode; 

float test; // soak up number 

long facevalue; 

long Unit, Atom; 

string Name; 

string Junk; 

long counter = 0 ; 

while (! ExpStream. eof () ) 

{ 

ExpStream >> facevalue >> TestString >> Name >> Junk >> Unit >> Atom; 

if (TestSt ring. length 0 >1 && Unit==UnitLimit ) 

{ 

NewNode = new ProbeNode; 
NewNode- >DataPointer = new ProbeSynth; 
NewNode- >DataPointer->Original = TestString; 
transform (TestString) ; 

NewNode->DataPointer->Probe = TestString; 
NewNode- >DataPointer->Name = Name; 
NewNode- >DataPointer->Locat ion = facevalue; 
NewNode- >DataPointer->Unit = Unit; 
NewNode- >DataPointer->Atora = Atom; 

NewNode- >DataPointer->SynthesizeMisMatch (MutateValue) ; 
this->DataSet . InsertAf terCurrent (NewNode) ; 

} 

counter ++; 



22 



if (counter%100==0) 
cout « counter << TAB << Name << TAB << TestString <:< endl; 
} 

ExpStream. close () ; 

cout << DataSet-Cost << endl; 

} 

TSPClass: :LoadExpressionFluff Data (string FileName, long MutateValue) 
{ 

// build the basic datalist 
ifstream ExpStream,- 

ExpStream. open (FileName .c_3tr (} , ios: : in) ; 

if (ExpStream.badO ) 

{ 

cout << endl << "Experimental data file specified does not exist!"; 
cout << endl << endl; 
exit (0) ; 

} 

long probeloop = 0; 

string TestString; 

ProbeNode *NewNode; 

float test; // soak up number 

long facevalue; 

long Unit, Atom; 

string Junk; 

string Name ; 

long counter = 0 ; 

while ( lExpStream.eof (J ) 

{ 

ExpStream » facevalue » TestString >> Name; 
if (TestString. length 0 >1) 
{ 

NewNode = new ProbeNode; 
NewNode->DataPointer = new ProbeSynth; 

transform (TestString) ; 
NewNode- >DataPointer->Probe = TestString; 
NewNode- >DataPointer->Name = Name; 

NewNode- >DataPointer->Locat ion = facevalue; 
NewNode->DataPointer->SynthesizeFastFluf f (MutateValue) ; 
this->DataSet - InsertAf terCurrent (NewNode) ; 

} 

counter++ ; 

if (counter%100==0) 

cout << counter << TAB << Name « TAB << TestString << endl; 

} 

ExpStream . close ( } ; 

cout << DataSet.Cost << endl; 

} 

TSPClass: :LoadChipData( string FileName) 
{ 

// build the basic datalist 
ifstream ExpStream; 

ExpStream. open (FileName .c_3tr 0 , ios: :in) ; 

if (ExpStream.badO) 
{ 

cout << endl << "Experimental data file specified does not exist!"; 
cout << endl << endl; 
exit (0) ; 

} 

long probeloop = 0; 
string TestString; 
ProbeNode *NewNode; 

long test; // soak up number 

while (! ExpStream. eof 0 ) 
{ 

ExpStream >> TestString >> test; 
if (TestString. length 0 >1) 
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{ 

NewNode = new ProbeNode; 
NewNode->DataPointer = new ProbeSynth; 

NewNode- >DataPointer->Probe = TestString; 
NewNode->DataPointer->SynthesizeMutant (test-l) ; 
this->DataSet . InsertAEterCurrent {NewNode) ; 

} 

} 

ExpStream. close {) ; 

cout « DataSet.Cost << endl; 

} 

TSPClass: -.OutputTour (string FileName) 
{ 

of stream OutStreatn; 

OutStream.open (FileName . c_str ( } , ios: :out) 
this->Be St Tour .Output (OutStream) ; 
OutStreatn. close 0 ; 

) 

TSPClass: :ListTour (string FileName) 
{ 

of St ream OutStream; 

OutStreatn. open (FileName . c_str () , ios: ;out) ; 
this->BestTour. List (OutStream) ; 
OutStream. close () ; 

} 

TSPClass: ;AppendTour (string FileName) 
{ 

ofstream OutStream; 

OutStream.open(FileName.c_str 0 , ios: :app) ; 
this->BestTour -Output (OutStream) ; 
OutStream- close { ) ; 

} 

static void do_gematch {long Value) 
{ 

TSPClass Example; 

Example. LoadMisMatchFluff Data ( "getnatch.prb" , 10-1) ; 
Example -GenerateTourBylnsertionWithModifier (Value) ; 
Exampl e . Outpu t Tour { "gemat ch . f I f " ) ; 

} 

static void do_yematch() 
( 

TSPClass Example; 

Example . LoadExpressionData ( "yematch . prb" , 10-1) ; 
Example .GenerateTourBy Insert ionWithModifier(l) ; 
Example . OutputTour ( "yematch . tsp" ) ; 

} 

static void do_yematch local (long UnitMatch) 
{ 

TSPClass Example; 

Example -LoadExpressionDataByUnit { "yematch . prb" , 10-1, UnitMatch) 
Example .GenerateTourBylnsertionWithModifier (1) ; 
Example -AppendTour ( "yematch . Isp" ) ; 

) 

static void do_yematch_local_pool (long UnitMatch) 
{ 
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TSPClass Example; 

Example. LoadExpressionDataByUnit {"yematch.prb" , 10-1, UnitMatch) ; 
Example. GenerateTourByClosest Pool () ; 
Example . AppendTour { "yematch - csp" ) ; 



static void do_hiimmatch_local {long UnitMatch) 
{ 

TSPClass Example; 

Example. LoadExpressionData ("ha. prb", 13-1) ; 
Example . Genera teQuickTourBy Insert ion (2048) ; 
Example . OutputTour ( "ha . tsp" ) ; 

//while (Example. GenerateTourBylnsertionAndDelet ion {) ) 
//Example . AppendTour { "ha . Isp" ) ; 

} 

static void do fast_gematch(long Value) 
{ 

TSPClass Example; 

Example. LoadMisMatchFastFluff Data ("gematch. prb" , 10-1) ; 
Example -GenerateTourBylnsertionWithModifier (Value) ; 
Exampl e. OutputTour ("gefast. flf ") ; 



static void do_f ast_pool_gematch(long Value) 
{ 

TSPClass Example; 

Example -LoadMisMatchFastFluff Data ( "gematch. tny" , 10-1) ; 
Example .GenerateTourByClosestPool () ; 
Example .OutputTour { "gefast .pi") ; 



static void do_mito(long Value) 
{ 

TSPClass Example; 

Example. LoadMutantData("mt9566 .prb", 10-1) ; 

Example .GenerateTourBylnsertionWithModifier (Value) ; 

/ /Example . ImproveTourByReplacement (1000) ; 

string Test; 

Test = "mt95S6."; 
char ctest = 'a' rvalue; 

Test += ctest; 
Example .OutputTour (Test) ; 

} 



static void do_hiv(long Value) 

{ 

TSPClass Example; 

Example . LoadChipData ( "hv43 0a .prb" ) ; 

Example .GenerateTourBylnsertionWithModifier (Value) ; 

//Example . ImproveTourByReplacement (1000) ; 

string Test ; 
Test = "hv430."; 

char ctest = 'a'+Value; 
Test += ctest; 

Example .OutputTour (Test) ; 

} 

static void do_new_ge() 
{ 

TSPClass Example; 
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Example .liOadExpressionFluff Data ( "i : \\data\\roetrix\\ehubbe\\design\\nmi3c\\gera01\\9eTn01 .prb", 

10-1) ; 

Example .GenerateTourBy Insert ionWithMod if ier (1) ; 
Example . OutputTour ( " geraO 1 . f 1 f " ) ; 

} 



static void do_noise () 
{ 

TSPClass Example; 

Example . LoadData { "noiseaS . 20 " ) ; 

Example . GenerateTourBylnsertionWithModif ier (1) ,- 

Example . ListTour ( "na8_20 . tsp" } ; 

} 

static void do_noise_two { ) 
{ 

TSPClass Example; 

Example . LoadData ( "c : \\cover\\noisea8 . 16" ) ; 
Example . GenerateTourBy Insert ionWithModif ier { 1 ) ; 
Example . ListTour ( "na8_16 . tsp" ) ; 

} 

static void do_noise_three { ) 

{ 

TSPClass Example; 

Example . LoadData ( "c : \\cover\\noisea8 . 18" ) ; 
Example -GenerateTourBylnsertionWithModif ier (1) ; 
Example . ListTour ( "na8_18 . tsp" ) ; 

} 

static void do_cost_one ( ) 
{ 

TSPClass Example; 

Example . LoadData ( "c : WgeniusWtestaB . 20") ; 
Example . GenerateTourBylnsert ionWithModif ier { 1 ) ; 
Example . ListTour ( "ca8_2 0 . tsp" ) ; 

} 

static void do_cost_two { } 
{ 

TSPClass Example; 

Example . LoadData ( "c : \\9enius\\ca8_20 . prb" ) ; 
Example .GenerateTourBylnsertionWithModif ier (1) ; 
Example -ListTour ( "cba8 20. tsp") ; 

} 

static void do_cost_quick { ) 
{ 

TSPClass Example; 

Example -LoadData ( "c : \\9enius\\noise\\ca8_20 . rnd" ) ; 
Example . GenerateQuickTourBylnsert ion (1024 ) ; 
Example .ListTour { "cqa8_20 .45" ) ; 

} 

static void do_rat_local (long UnitMatch) 
{ 

TSPClass Example; 

Example .LoadExpressionData ( "r : \\alldes\\cdesign\\ter09\\included_probes\\normal_probes . txt " , 
13 -1 ) ; 

Example -LoadExpressionData {"r: \\alldes\\cdesign\\ter09\\included_probes\\sense\\rever3e_comp 
_seqs.txt" , 13 -1) ; 

// Example .LoadSingleExpressionData {"r: \\alldes\\cde3ign\\eol0191\\eoshu02 .dat ") ; 

// Example -GenerateQuickTourBylnsertion {140000) ; 
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Example .GenerateQuickTourByInsertion{20480) ; 

Example .OutputTour { "r : \\alldes\\cdesign\\ter09\\full_£ix . tsp" ) , 

//while (Example.GenerateTourBylnsertionAndDeletionO ) 
//Example . AppendTour ( "ha . Isp" ) ; 



main ( ) 
{ 

do_rat_local { 0 ) ; 

// do_f ast_pool_gematch(l) ; 

// do_yeniatch() ; 

// for (long i=l; i<104 ; i++) 

// do_yematch_local_pool (i) ; 

// do_gematch(l) ; 

// do_new_ge 0 ; 

// do_noise ( ) ; 

// do_noise_two ( ) ; 

// do_noise_three ( ) ; 

// do_cost_quick 0 ; 



} 
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// Edge Optimizer 
#include <cstring.h> 
#include <f stream. h> 
#include <st(iio.h> 
#include <stdlib.h> 
#include <dir.h> 
#include <time.h> 
#include <math.h> 

#define TAB ' \t ' 
#define MXNAME 40 
#define MXLINE 1000 

#define TRUE 1 
#define FALSE 0 
#define MXSEQ 4 5 
#define MXFEATURE 4 5 
ftdefine MXQUALIFIER 45 

// Edge Optimizer Story 

// 1) Strip off all Valid Blocks 

// 2) Put Valid Blocks On to Minimize Edges 

// Input: Cdl file. Ret file. Parameters 
// Output : Twisted Cdl file and Ret file 

char complement (char base) 
{ 

if {base=='A') 

return ( *T' ) ; 
if (base=='C') 

return TG ' ) ; 
if (base=='G') 

return ( ' C ' ) ; 
if {base=='T') 

return ( 'A' ) ; 
return (base) ; 

} 

class EntryClass{ 

// what CDL information is associated with everything 
public : 

char sequence [MXSEQ] ; 
int destype; 

char feature tMXFEATURE] ; 



char qualifier [MXQUALIFIER] ; 

int expos; 

int endpos; 

int pos; 

char pbase [MXFEATURE] , tbase [MXFEATURE] ; 
int finishpos; 
int fixed; 
int variable; 

int unit, block; 
long atom; 
int repeat ; 
int seqno ; 
long layout; 

char locus [MXFEATURE] ; 
char accession [MXFEATURE] ; 
EntryClassO { Initialize ( ) ; } ; 
Initialize () ; 
LineScan (char *) ; 
DumpLine (FILE *fp, int i, int j); 
DumpHut (FILE *£p) ; 



EntryClass : : Initialize () 
{ 

strcpy (sequence, ""); 

destype = 0; 

strcpy ( feature, ""); 

strcpy (qualifier , ""); 
expos = 0; 
pos = 0; 

strcpy (pbase ,"!"); 
strcpy (tbase, " ! " ) ; 
unit = 0 ; 
block = 0; 
atom - 0; 

} 

EntryClass :: LineScan (char *Line) 
{ 

int X,Y; 

static char PROBE [MXLINE] ; 
int DESTYPE; 

Static char FEATURE [MXLINE] , 

QUALIFIER [MXLINE] ; 

int EXPOS ; 

char TBASE [MXLINE] ; 

int ENDPOS, 

POSITION; 
char PBASE [MXLINE] ; 
int FINISHPOS, 

FIXED, 

VARIABLE , 

UNIT, 

BLOCK, 

REPEAT, 

SEQNO, 

LAYOUT; 

long ATOM; 

static char ACCESSION [MXLINE] , 
LOCUS [MXLINE] ; 

sscanf(Line, "%d %d %s %d %s %3 %d %s %d %d %s %d %d %d %d %d %ld %d %d %d %s 

%s", 

&X, 
&Y, 

PROBE , 

&DESTYPE, 

FEATURE, 

QUALIFIER, 

tEXPOS, 



TBASE, 

&ENDPOS, 

&POSITION, 

PBASE , 

&FINISHPOS, 

&FIXED, 

^VARIABLE, 

&UNIT, 

&BLOCK, 

&ATOM, 

&REPEAT, 

&SEQNO, 

&LAYOUT, 

LOCUS, 

ACCESSION) ; 

St rcpy (sequence, PROBE); 

destype = DESTYPE; 

strcpy (feature, FEATURE) ; 

strcpy(qualif ier, QUALIFIER) ; 

strcpy (locus, LOCUS) ; 

strcpy (accession, ACCESSION) 

expos = EXPOS ; 

endpos = ENDPOS; 

finishpos = FINISHPOS; 
fixed = FIXED; 
variable = VARIABLE; 

repeat = REPEAT; 

segno = SEQNO; 

layout = LAYOUT; 

Strcpy (tbase, TBASE); 
strcpy (phase, PBASE) ; 
unit = UNIT; 
block = BLOCK; 
atom = ATOM; 
pos = POSITION; 

) 

EntryClass : rDutnpLine (FILE *fp, int i, int j) 
{ 

fprintf (fp, 

"%d\t%d\t%s\t%d\t%s\t%s\t%d\t%s\t%d\t%d\t%s\t%d\t%d\t%d\t%d\t%d\t%ld\t%d\t%d\t%d\t%s\t%s\: 

i , 
j, 

sequence, 

destype , 

feature, 

qualifier, 

expos , 

tbase , 

endpos , 

pos, 

phase, 

finishpos, 

fixed, 

variable, 

unit, 

block, 

atom, 

repeat , 

seqno, 

layout, 

locus, 

accession) ; 

} 

EntryClass : rDumpMut (FILE *fp) 
{ 

char tempc; 



if (destype==0) 



{ 

fprintf(fp, "-"); 
return (TRUE) ; 

) 

if {abs (deatype) <:100) 
{ 

if (destypoO) 

fprintf(fp, "C"}; 

else 

fprintf(fp, "X"); 
ret urn (TRUE) ; 

} 

if (strlen(pbase) >1) 

{ 

f print f(fp, "I"); 
return (TRUE) ; 

} 

if (phase [0]==' ! •) 
{ 

Eprintf{fp, "D"); 
ret:urn(TRUE) ; 

} 

if fdestype>0) 
{ 

tempc = complement (pbase [0] ) ; 

} 

else 

tempc = pbase [0] / 
if (tempc==tbase [0] ) 
{ 

fprintf(fp, "%c", tempc); 
return (TRUE) ; 

} 

fprintf (fp, " ") ; 
return (TRUE) ; 

} 



class SynthClass( 

// how things are built 

public : 

char * synthesis; 

int synlength; 

SynthClass () {synthesis=NULL; synlength = 0;}; 

Allocate (int) ; 

DeAllocate 0 ; 

Dif f (SynthClass &) ; 

SetBit (char, int) ; 

char Get Bit (int) ; 

GetLast ( ) ; 

GetFirst (} ; 

-SynthClass 0 ; 



SynthClass: : -SynthClass () 
{ 

DeAllocate 0 ; 

} 

SynthClass :: Allocate (int Size) 
{ 

synthesis = new char [Size] ; // just a bitfield, really 
if (synthesis==NULL) 
{ 

printfC'Blow upi In SynthClass :: Allocate") ; 
return (FALSE) ; 

} 

for (int i=0; i<Size; i++) 

synthesis [i] =0; // Nothing, null, no bits set! 
synlength = Size; 

return (TRUE) ; 




} 

SynthClass: :DeAllocate {) 
{ 

if (synthesis ! =NULL) 

deleted synthesis; 
synlength = 0; 
synthesis = NULL; 

} 

SynthClass :: SetBit (char value, int Which) 
{ 

if (synthesis 1 =NULL Which< synlength && Which>-1) 
synthesis (Which] = value; 

} 

char 

SynthClass: ;GetBit (int Which) 
{ 

if (synthesis! =KULL && Which< synlength && Which>-1) 
return (synthesis {Which] ) ; 

else 

return (0) ; 

1 

SynthClass: :Diff (SynthClass StSource) 

^0 { 

ifi int count = 0; 

if (! (synthesis ! =NULL Source . synthesis ! =NULL synlength==Source . synlength) } 
■?=* return (0) ; 

O for (int i=0; i<synlength; i++) 

count += (synthesis [i] ! =Source . synthesis [i] ) ; 
: = return ( count ) ; 

T ^ 

'3 SynthClass i : Get Last ( ) 

m ^ 

for (int i=synlength-l; i>0; i--} 
i= if (synthesis [i] >0) 

ij\ return (i) / 

= 2 SynthClass : : GetFirst ( ) 
{ 

for (int i=0; ioynlength; i++) 
if (synthesis [i] >0} 

return (i) ; 

} 

class LocalDataClass{ 
public : 

int relx, rely; // relative positions within a block 
int validflag; // newly added - can I move this for optimization? 

EntryClass cdldata; 
SynthClass retdata; 

LocalDataClass 0 {relx=rely=0; validflag = TRUE;); 
-LocalDataClass 0 ; 
SetRelative (int , int); 
PrintLcngSeq (FILE *fp) ; 

}; 

LocalDataClass : ; -LocalDataClass ( ) 
{ 

retdata. DeAllocate {) ; 

) 

LocalDataClass :: SetRelative (int i, int j) 
{ 

relx = i; 



rely = j; 

) 

LocalDataClass : : Print LongSeg (FILE *fp) 
{ 

int j, i=3; 

for (j=0; j<retdata. synlength; j++) 

{ 

if (retdata.GetBit (j) ) 
{ 

fprintf(fp, "%c", cdldata. sequence ti] ) ; 
i++; 

} 

else 

fprintf(fp, "."); 

) 

) 



class BlockClass{ 
public : 

LocalDataClass **DataStack ; 
int DataStackSize ; 

int WSize, HSize; // rectangular grid 

BlockClass 0 {DataStack=NULL; DataStackSize = 0; WSize = HSize = 0;} 
Allocate (int) ; 
DeAl locate {) ; 
-BlockClass () ; 

); 

BlockClass :: Allocate (int Size) 
{ 

DataStack = new LocalDataClass * [Size] ; 
if (DataStack==NULL) 

return (FALSE) ; 
for {int i=0; i<Size; i++) 

DataStack [i] =NULL; 
DataStackSize = Size; 



BlockClass: ;DeAllocate () 
{ 

if (DataStack==NULL) 
return (TRUE) ; 
for (int i=0; i<DataStackSize; i++) 
if (DataStack [i] !=NULL) 

{ 

printf ( "Error ! Data Leakage from Block ! \n" ) ; 
delete DataStack til ; 
DataStack [i] = NULL; 

} 

delete n DataStack; 
DataStack = NULL; 
DataStackSize = 0; 
WSize = HSize = 0; 



BlockClass: : -BlockClass {) 
{ 

DeAllocate 0 ; 

} 

class BlockStackClass{ 
public: 

BlockClass **BlockStack; 
long BlockCurSize; 
long BlockStackSize; 

BlockStackClass () {BlockStack=WLL; BlockStackSize = 0 ; BlockCurSize 
Allocate (long) ; 



DeAllocateO ; 

-BlockStackClassO (DeAllocateO } ; 
BlockClass *PutBlockOnStack (BlockClass *) ; 
BlockClass * Remove Block (long) ; 
BlockClass *TeraporaryBlockFromStack ( long) 
Swap {long, long) ; 
Shuffle 0 ; 

}; 

BlockStackClass : : Allocate (long Size) 
{ 

BlockClass ** TmpStack; 

TmpStack = new BlockClass * [Size} ; 
if (TnipStack==NULL) 

return (FALSE) ; 
long i; 



for {i=0; i<Size; i++) 
TtnpStack [il =NULL; 

for (i=0; i<BlockCurSi2e i<BlockStackSize && BlockStack!=NULL; i++) 
{ 

TmpStackti] = BlockStack [i] ; 
BlockStackti] = NULL; 

} 

if (BlockStack !=NULL) 

delete [] BlockStack; 
BlockStack = TmpStack; 
TmpStack = NULL; 
BlockStackSize = Size; 
return (TRUE) ; 



BlockClass * 

BlockStackClass: : PutBlockOnStack (BlockClass *TempBlock) 



if {BlockCurSize<BlockStackSize) 

{ 
} 

else 
{ 

if { ! Allocate (BlockStackSize+1000) ) 

{ 

print f ("Can' t increase block stack\n") 
return (TempBlock) ; // upgrade stack 

} 

} 

BlockStack [BlockCurSize] = TempBlock; 
BlockCurSize++ ; 

return (NULL) ; // remove pointer 



BlockStackClass :: Swap (long Source, long Sink) 
{ 

BlockClass *TeTnpBlock; 



TempBlock = BlockStack [Sink] ; 
BlockStack [Sink] = BlockStack [Source] ; 
BlockStack [Source] = TempBlock; 
TempBlock = NULL; 
return (TRUE) ; 



BlockClass * 

BlockStackClass: :RemoveBlock(long whichBlock) 
{ 

BlockClass *TempBlock; 



if (WhichBlock>=BlockCurSi2e | | WhichBlock<0) 
return (NULL) ; 



TempBlock = BlockStack IWhichBlock] ; 
BlockCurSize-- ; 

BlockStack [whichBlock] = BlockStack [BlockCurSize] ; 
BlockStack [BlockCurSize] =NULL; 
return (TempBlock) ; 

} 

BlockClass * 

BlockStackClass : :TemporaryBlockFromStack{long WhichBlock) 
{ 

BlockClass *TempBlock; 

if (WhichBlock>=BlockCurSize | | WhichBlock<0) 
return (NULL) ; 

TempBlock = BlockStack [WhichBlock] ; 
return (TempBlock) ; 

} 

BlockStackClass: iDeAllocate ( ) 
{ 

if (BlockStack==NULL) 
return (TRUE) ; 
long i; 

for {i=0; i<BlockStackSize; i++) 
{ 

if (BlockStack [i] ! =NULL) 
{ 

printfC'Data Leakage from BlockStack\n" ) ; 
delete BlockStack [i] ; 
BlockStack [i] = NULL; 

} 

} 

delete [] BlcckStack; 
BlockStack = NULL; 
BlockStackSize = 0; 
BlockCurSize = 0; 

} 

BlockStackClass: : Shuffle () 
{ 

long t ; 
long i; 

// randomly rearrange the stack to prevent bias 
for (i=BlockCurSize-l ; i>0; i--) 
{ 

t = randO %32000; 
t = t*32000+(rand{) %32000) ; 
t = t % (i+1) ; 

Swapd, t) ; 

} 

} 

class ChipArrayClass{ 
public : 

BlockStackClass ValidBlockStack; 
LocalDataClass ***DataGrid; 
int Xdim; 
int Ydim; 
int SynthSteps; 

char ret name [MXNAMEJ ; 

long NumOnes; // useful statistic 

int GlobalHeight ; 
int MaxAllowed; 
int Radius ; 



float LeakageHalfLife; 
int ScanRadius; 



int weightflag; // type of weights to use 



// constraints 



ChipArrayClass () ; 
-ChipArrayClass ( ) ; 
Allocate {int , int) ; 
DeAllocate ( ) ; 
ReadcdKchar *, int) ; 
DunipCdl (char *) ; 
ReadRet(char *, int) ; 
DumpRet {char *) ; 

StripAreaToBlock (BlockClass * , int , int , int , int) ; 
CheckBloclcFitToArea (BlockClass *, int, int) ; 
PutBlocklnArea (BlockClass *, int, int) ; 

ValidMove (int , int) ; 
ValidLocation (int , int) ; 
Valid (int, int) ; 
Val idBlock { int , in t , int , int ) ; 
ValidTile { int , int , int , int ) ; 
ValidBlank (int , int, int, int) ; 

=1 CountDif f (LocalDataClass *, int, int); 

g double CountEdges {BlockClass *, int, int) ; 

double CountWeightedEdges (BlockClass *, int, int) ; 

double Count EdgesFromStack (long, int, int) ; 

FindNextDiagonalSlot (int &, int &) ,- 
== FindNextHorizontalSlot (int &, int &) ; 

^ DiagonalReplacement (long) ; 

^ HorizontalReplacement (long) ; 

Ti StripAllValidBlocks (int) ; 

=; PlaceBlockPromStack (long, int, int), - 

- SearchLocationWithStats (int , int , long, long &, double &, double &, double &) 

2 Count UnitlnArea {long, int, int , int , int ) ; 

^ ProximityCheckBlock (BlockClass *, int , int , int , int) ; 

^ ProximityCheckFromStack ( long , int , int , int , int ) ; 

StripBadProximityValues (int ) ; 
1; PickRandomValidBlock ( int &, int int, int) ; 

^ ReadlnstructionFile (char *); 

~ InterpretlnstructionLine (char *) ; 

GenerateMutFile (char *) ; 

SetUnits (long, long, int) ; 

SetArea(int, int, int, int, int); 

SetAnt iArea ( int , int , int , int , int ) ; 

SetDestype (int , int); 

Shuffle 0 ; 

StripValidBlock(int, int, int) ; 
StripRandomBlocks (long, int); 
DoubleReplacement (long) ; 
GenerateDif EFile (char *) ; 



FindNextAggregateSlot (int &, int &) ; 
AggregateReplacement (long) ; 



ChipArrayClass: :Shuffle() 
{ 

ValidBlockStack.ShuffleO ; 

} 

ChipArrayClass :: Set Area (int X, int Y, int tX, int tY, int value) 
{ 

int i , j ; 

for {i=X; i<=tX; i++) 



for (j=Y; j<=tY; 
{ 

if {valid(i,j)) 

DataGrid [i] [ j] ->validf lag = value; 

} 

} 

} 

ChipArrayClass : : SetAntiArea (int X, int Y, int tX, int tY, int value) 
{ 

// used to set all >but< a given physical area 
int i , j ; 

for (i=0; i<Xditn; i++) 
{ 

for (j = 0; j<YdiTn; + ) 
{ 

if (Valid(i,j)) 

{ 

if (! (i>=X && i<=tX j>=Y && j<=tY)) 

DataGrid[i] [ j ] - >validf lag = value; 

} 

} 

} 

} 

ChipArrayClass :: SetDestype {int destype, int value) 
{ 

int i,j; 

for (i=0; i<Xdim; i++) 



{ 



for (j=0; j<Ydim; 
{ 

if (Valid (i,j)) 

{ 

if (destype==DataGrid [i] [j ] ->cdldata . destype) 
{ 

DataGrid[i] [ j ] - >validf lag = value; // can't move this one 

} 



ChipArrayClass: :ValidMove{ int X, int Y) 
{ 

if (DataGrid [X] [Y] ->validf lag) 
return (TRUE) ; 
return (FALSE) ; 



ChipArrayClass : :ValidLocat ion (int X, int Y) 
{ 

if (X<0 I I Y<0 1 I X>=Xdim | j Y>=Ydim) 
return (FALSE) ; 
// nothing here yet 
return (TRUE) ; 

} 

ChipArrayClass :: Valid (int X, int Y) 
{ 

if { !ValidLocation(X,Y) ) 

return (FALSE) ; // not a location we're allowed to move 
if (DataGrid [X] [Yl==mJLL) 

return (FALSE) ; // doesn't even exist! 
return (TRUE) ; 

} 

ChipArrayClass : :SetUnit3 (long Start, long Finish, int value) 
{ 



int i , j ; 

£or (i=0; i<Xdim; i++) 

for (j=0; j<Ydim; j++) 
{ 

if {Valid{i,j)) 

{ 

if (Start<=DataGrid [i] [j] ->cdldata.vinit && Finish>=DataGrid [i] [j ] ->cdldata.unit) 
{ 

DataGridfi] [j J ->validf lag = value; // can't move this one 



ChipArrayClass: :ValidTile(int X, int Y, int Width, int Height) 
{ 

long U, A; 
int tx, ty, i, j ; 
int xl, yl; 
xl = X+Width; 
yl = Y+Height; 

U = DataGrid[X] [Y] ->cdldata .unit ; 
A = DataGridLX] [Y] - >cdldata . atom; 
// if (Width==Height Height==l) // if we're moving control probes around! 
// return (TRUE); 

for (i=-l; i<width+l; i++) 

{ 

for (j=-l; j<Height+l; j++) 
{ 

tx = i+X; 

ty = j+Y; 

if {tx>=X ty>=Y && tx<xl && ty<yl) 
{ 

if ( ! valid(tx, ty) || ! ValidMove ( tx, ty) ) // if not valid we're unhappy 
return (FALSE) ; 

if (U!=DataGrid[tx] tty] ->cdldata.unit || A! =DataGrid [tx] [ty] - >cdldata . atom) // if 
not same , we ' re unhappy 

return (FALSE) ; 

} 

else 
{ 

if (Valid(tx, ty) ) 
{ 

if (U==DataGrid [tx] [tyj ->cdldata .unit && A==DataGrid ftx] [ty] ->cdldata . atom) 
// if outside same block unhappy 

return (FALSE) ; 

} 

) 

} 

} 

return (TRUE) ; 

} 

ChipArrayClass: :ValidBlank( int X, int Y, int Width, int Height) 
{ 

int i,j, tx,ty; 

// valid set of blanks 
for (i=0; i<Width; i++) 

for {j=0; 3<Height; j++) 
{ 

tx = i+X; 
ty = j+Y; 

if (! Valid (tx,ty) || ! ValidMove (tx, ty) ) 

return (FALSE) ; 
if {DataGrid[tx] [tyj ->cdldata .destype ! =0) 



return (FALSE) ; 

} 

return (TRUE) ; 

} 

ChipArrayClass: :ValidBlock(int X, int Y, int Width, int Height) 
{ 

if ( !Valid{X,Y) ) 

return (FALSE) ; 
if {]ValidMove(X,Y)) 

return (FALSE) ; 
if (ValidTile(X,Y, Width, Height)) 

return (TRUE) ; 

if (ValidBlank(X, Y, Width, Height)) // allow blank blocks to be moved, if validflag set for 
destype 0 

return (TRUE) ; 

return (FALSE) ; 

} 

ChipArrayClass :: PickRandomValidBlock( int &X, int &Y, int Width, int Height) 
{ 

long counter = 0 ; 
long Limit = 10000; 

X = rand()%Xdim; 
Y = randO %Ydim; 

while { !ValidBlock(X,Y, Width, Height) && counter< Limit) 

{ 

X = randO %Xdim; 
Y = randO %Ydim; 
count er++; 

} 

if Ccounter==Limit) 

return (FALSE) ; // can't find one in reasonable time! 
return (TRUE) ; 



ChipArrayClass : :CountDiff (LocalDataClass *TestData, int X, int Y) 
{ 

if ( !Valid(X, Y) ) 

return(O); // doesn't exist or is off chip, so no problems! 
return (TestData- >retdata. Dif f (DataGrid [X] [Y] ->retdata) ) ; 

} 

double 

ChipArrayClass : :CountEdges (BlockClass *TempBlock, int X, int Y) 
{ 

int i,tx,ty; 
int count = 0 ; 

// count the edges, if this block is in this location 
// note that "interior" edges of blocks are >not< counted 

for (i=0; i<TempBlock->DataStackSize; i++) 

{ 

if (TempBlock->DataStack[i] !=NULL) 
{ 

tx = X+TempBlock->DataStack [i] ->relx; 
ty = Y+TempBlock->DataStack [i] ->rely ; 
count +=CountDif f (TempBlock- >DataStack [i] , tx+1, ty) 
count +=CountDiff (TempBlock- >DataStack [i] , tx-i, ty) 
count +=CountDiff (TempBlock- >DataStack [i] , tx, ty-1) 
count +=CountDiff (TempBlock- >DataStack [i] , tx, ty+1) 

) 

} 

return (count ) ; 



doiible 

ChipArrayClass : iCountWeightedEdges (BlockClass *TempBlock, int X, int Y) 
{ 



int i,tx,ty; 
int rangex, rangey; 
doiible count = 0 ,- 
double Icount ; 
double distance; 

// count the edges, if this block is in this location 
// note that "interior" edges of blocks are >not<: counted 

for (i=0; i<TempBlock->DataStackSize; i++) 

{ 

if {TempBlock->DataStack[i] !=NULL) 
{ 

for (rangex = -l*ScanRadiu3 ; rangex<=ScanRadius; rangex++) 



{ 



for (rangey=-l*ScanRadius; rangey<=ScanRadius; rangey++) 
{ 

if (rangex!=0 (( rangey.'=0) 

{ 

distance = { rangex* rangex) + (rangey* rangey ) ; 
distance = sqrt (distance) ; 

tx = X+TempBlock- >DataStack [i] ->relx+rangex; 
ty = Y+TempBlock->DataStack[i] ->rely+rangey; 
Icount =CountDif f (TempBlock->DataStack [i] , tx, ty) 
Icount* =pow ( . 5, LeakageHalfLife*distance) ; 
count +=lcount; 



return (count) ; 



ChipArrayClass: : Count Unit InArea (long Unit, int X, int Y, int Width, int Height) 

{ 

int i , j ; 
int tx, ty; 
int count =0; 

long tatom=-100; // unlikely in any real unit 
for (i=0; i<Width; i++) 

{ 

tatom = -100; // start over with each vertical stripe - assumes vertical "atoms" 
for (j=0; j<Height; j++) 
{ 

tx = X+i; 
ty= Y+j; 

if (Valid{tx, ty) ) 
{ 

if (DataGrid[tx] [ty] ->cdldata.unit==Unit) 
{ 

// works because we're scanning vertically 

if (tatom! =DataGr id [tx] [ty] ->cdldata . atom) // only count a given unit/atom 

once 

count ++ ; 

tatom = DataGrid[txl [tyj ->cdldata . atom; 



} 

return (count) ; 

} 

ChipArrayClass: :ProxiTnityCheckBlock(BlockClass *TempBlock, int X, int Y, int Width, int Height) 
{ 

return (Count Unit InArea (TempBlock- >DataStack [0] - >cdldata . unit , X, Y, Width, Height ) ) ; 

} 

ChipArrayClass: :ProximityCheckFromStack( long Which, int X, int Y, int Width, int Height) 
{ 

BlockClass *TempBlock; 



int count; 

TempBlock = ValidBlockStack .TeraporaryBlockFromStack (Which) ; 
if (TempBlock==NULL) 

return (16000) ; // big error 
count = ProximityCheckBlock (TempBlock, X, Y, Width, Height) ; 
TempBlock = NULL; 
return (count) ; 

} 

double 

ChipArrayClass :: Count EdgesFromStack (long Which, int X, int Y) 
{ 

BlockClass *TempBlock; 
double count; 

TempBlock = ValidBlockStack . TemporaryBlockFromS tack (Which) ; 
if (TempBlock==NULL) 

return (16000) ; // big error 
if { iweightf lag) 

count = Count Edges (TempBlock, X, Y) ; 

else 

count = CountWeightedEdges (TempBlock, X,Y) ; 
TempBlock = NULL; 
return (count ) ; 

} 

ChipArrayClass :: FindNextDiagonalSlot (int &X, int &Y) 

{ 

// look for a NULL entry 
while (Y<Ydim && DataGrid [X] [Y] !=NULL) 
{ 

X-=l; 
Y+ = l; 

if (X<0 I j Y>=Ydim) 
{ 

X = Y+X+1; 
Y = 0; 

} 

if (X>=Xdim) 
{ 

Y = X- (Xdim-1) ; 
X = Xdim-1; 

} 

} 

if {Y==Ydim) // ran out of room! 

return (FALSE) ; 
else 

return (TRUE) ; 

} 

ChipArrayClass :: FindNextHorizontalSlot (int &X, int &Y) 
{ 

// look for a NULL entry 
while (Y<Ydim DataGrid [X] [Y] ! =NULL) 
{ 

X +=1; 

if (X>=Xdim) 
{ 

X = 0; 
Y++; 

} 

} 

if (Y==Ydim) // ran out of room! 

return (FALSE) ; 
else 

return (TRUE) ; 

} 

ChipArrayClass: : FindNextAggregateSlot (int &X, int &Y) 
{ 



} 



int total; 

total = X; 
if (Y>total) 

total = Y; 
while (Y<Yditn && DataGrid[Xl [Y] ! =NULL) 
{ 

// look for slots 
if (X>Y) 

Y++; // move vertically 

else 

if (X<=Y) 

X--; // basic zero moves 

if (X<0) 
{ 

X=Y+1; // add one to total 
Y=0; 

} 

if (X>=Xdim) 
{ 

Y = X; 
X=Xdim-l; 



if (Y>=Ydim) 

return (FALSE) ; 
else 

return (TRUE) ; 



ChipArrayClass : : PlaceBlockFromStack (long Which, int X, int Y) 
{ 

BlockClass *TempBlock; 

TerapBlock = ValidBlockStack . RemoveBlock (Which) ; 
if (TempBlock J =NULL) 
{ 

if (PutBlocklnArea (TempBlock, X,Y)) 
TempBlock = NULL; // keep wacky pointers from drifting 

else 
{ 

printf ("Failure to fit block in area: %d %d\n", X, Y) ; 
TempBlock = ValidBlockStack . PutBlockOnStack (TempBlock) ; // throw back on stack 

} 



} 



else 

printf ("Failure to get from stackj %d %d\n", X,Y),- 



ChipArrayClass : : SearchLocationWithStats (int X, int Y, long searchlimit, long &Best, double &bestc, 

double &avg, double &worstc) 

{ 

long search; 
long count = 0 ; 
doiible c; 

Best = ValidBlockStack .BlockCurSize-i; // which one 

bestc = CountEdgesFromStack{ValidBlockStack.BlockCurSize-l,X, Y) ; 

avg = bestc; 
worstc = bestc; 
count = 1; 

for {searchsl; searchoearchlimit && 3earch<ValidBlockStack.BlockCurSi2e; search++) 
{ 

c = CountEdgesFromStack { ValidBlockStack. BlockCurSize-1- search, X,Y) ; // what 

if we put here? 

avg +=c,- 

if (c<:bestc) 

{ 

bestc - c; 

Best = ValidBlockStack. BlockCurSizG-l- search; 

} 

if (oworstc) 



worstc=C; 

} 

count ++ ; 

} 

avg /=count; // number actually searched 
return (TRUE) ; 



ChipArrayClass: : StripBadProximityValues {int H) 
{ 

int i , j ; 
long U; 
int c; 

BlockClasa *TerapBlock; 
GlobalHeight = H; 

for (i=0; i<Xdim; i++) 

{ 

for {j=0; j<Ydim; 

{ 

if (ValidBlock(i, j , 1, H) ) // note that blanks could mess this up badly! 

{ 

U = DataGrid[i] [ j ] ->cdldata .unit ; 

c = CountUnitInArea{U, i -Radius, j-Radius, 2*Radius+l, 2*Radius+l) ; 
if (oMaxAl lowed) 

StripValidBlock{i, j ,H) ; 

} 

} 

} 

// now we've got all our trovible removed from the chip 
printfC'Bad Proximity values: %d %d %d %ld\n", H, Radius, MaJxAllowed, 
ValidBlockStack.BlockCurSize) ; 
return (TRUE) ; 

} 



ChipArrayClass: : DoubleReplacement (long searchlimit) 
{ 

// idea is to "dilute" any bad values 
// this only works if the chip is sufficiently large 
// and there are sufficiently few bad items 

StripBadProximityValues (GlobalHeight) ; 

while (ValidBlockStack . BlockCurSize>0) 

{ 

StripRandomBlocks (ValidBlockStack.BlockCurSize+100, GlobalHeight); // get some good random 
locations freed up 

Shuffle (}; // rearrange life 

DiagonalReplaceraent (searchlimit) ; // put 'em back, - if too much search, goes back exactly 
to bad spots 

StripBadProximityValues (GlobalHeight) ; // find out if we've got them all 

) 

) 



ChipArrayClass: : DiagonalReplaceraent (long searchlimit) 
{ 

// replaces blocks from stack onto the chip 
int X, Y; 
long Best; 

double EdgesAdded = 0.1; 
double TotalAdded = 0.1; 
double AvgEdges =0.1; 
double WorstEdges=0 . 1 ; 
double bestc; 
dovible avg; 



do\able worstc; 

long Report = 1000000; 

Report /=3earchlimit ; 
if (Report>1000) 
Report =1000; 

X=Y=0; 

while {FindNextDiagonalSlot (X,Y) ) 
{ 

// found a location where a block was removed 

SearchLocationWithStats (X, Y, searchlimit , Best, bestc, avg, worstc); 
// found the best thing to put there 

// and so put it there ! 
PlaceBlockFromStackCBest,X,Y) ; 

EdgesAdded += bestc; AvgEdges += avg; WorstEdges += worstc; TotalAdded ++; 

if (ValidBlockStack . BlockCurSize%Report==0) 

printfC'At: %d %d %ld %lf %lf %lf\r", X, Y, ValidBlockStack. BlockCurSize 
EdgesAdded/ (2 *TotalAdded) , EdgesAdded/AvgEdges, EdgesAdded/WorstEdges) ; 
} 

printf ("\nAt : %d %d %ld %lf %lf %lf\n", X, Y, ValidBlockStack . BlockCurSize, 
EdgesAdded/ (2*TotalAdded) , EdgesAdded/AvgEdges, EdgesAdded/WorstEdges} ; 
return (TRUE) ; 

} 

ChipArrayClass : : AggregateReplacement (long searchlimit) 
{ 

// replaces blocks from stack onto the chip 
int X, Y; 
long Best; 

double EdgesAdded = 0 . l ; 
double TotalAdded = 0.1; 
double AvgEdges =0.1; 
double WorstEdges=0 . 1 ; 
double bestc; 
double avg; 
double worstc; 
long Report = 1000000; 

Report /=searchlimit ; 
if (Report>1000) 
Report =10 00 ; 

X=Y=0; 

while (FindNextAggregateSlot (X, Y) (ValidBlockStack . BlockCurSi2e>0) ) 
{ 

// found a location where a block was removed 

SearchLocationWithStats {X, Y, searchlimit , Best, bestc, avg, worstc); 
// found the best thing to put there 

// and so put it there! 
PlaceBlockFromStack(Best ,X, Y) ; 

EdgesAdded += bestc; AvgEdges += avg; WorstEdges += worstc ; TotalAdded ++; 

if (ValidBlockStack . BlockCurSize%Report==0) 

printfC'At: %d %d %ld %lf %lf %lf\r", X, Y, ValidBlockStack. BlockCurSize 
EdgesAdded/ (2*TotalAdded) , EdgesAdded/AvgEdges, EdgesAdded/WorstEdges) ; 
} 

printf ("\nAt: %d %d %ld %lf %lf %lf\n", X, Y, ValidBlockStack .BlockCurSize, 
EdgesAdded/ (2*TotalAdded) , EdgesAdded/AvgEdges, EdgesAdded/WorstEdges) ; 
return (TRUE) ; 

} 



ChipArrayClass: rHorizontalReplacement (long searchlimit) 



// replaces blocks from stack onto the chip 
int X, Y; 
long Best; 

dovible EdgesAdded = 0.1; 
double TotalAdded = 0.1; 
double AvgEdges =0.1; 
double WorstEdges=0. 1; 
doiible bestc; 
double avg; 
double worstc; 
long Report = 1000000; 

Report /ssearchlimit ; 
if (Report?! 000) 
Report =10 00; 

X=Y=0; 

while (FindNextHorizontalSlot (X, Y) ) 
{ 

// found a location where a block was removed 

SearchLocationWithStats(X, Y, searchlimit, Best, bestc, avg, worstc); 
// found the best thing to put there 

// and so put it there! 
PlaceBlockFroTtiStack(Best,X, Y) ; 

EdgesAdded += bestc; AvgEdges += avg; Worst Edges += worstc; TotalAdded++ ; 

if {ValidBlockStack.BlockCurSize%Report==0) 

printfC'At: %d %d %ld %lf %lf %lf \r" , X, Y, ValidBlockStack . BlockCurS 
EdgesAdded/ (2*TotalAdded) , EdgesAdded/AvgEdges, EdgesAdded/worst Edges) ; 
} 

printf ("\nAt: %d %d %ld %lf %lf %lf\n", X, Y, ValidBlockStack .BlockCurSize, 
EdgesAdded/ (2 *TotalAdded) , EdgesAdded/AvgEdges , EdgesAdded/WorstEdges) ; 
return (TRUE) ; 

} 

ChipArrayClass : : StripValidBlock (int X, int Y, int H) 
{ 

BlockClass *TeTnpBlock; 

if (ValidBlock(X,Y, 1, H) ) 
{ 

TempBlock = new BlockClass; 

St ripAreaToB lock {Temps lock, X,Y,1,H); // take probe from chip 
TempBlock = ValidBlockStack . PutBlockOnStack (TempBlock) ; 
if (TempBlock !=NULL) 
{ 

printf ("Failure to put on stackl %d %d\n" , X,Y} ; 

} 

} 

) 

ChipArrayClass : : StripAllValidBlocks ( int H) 
{ 

int i,j; 
BlockClass *TempBlock; 

for (i=0,- i<Xdim; i++) 

{ 

for (j=0; j<Ydim; 
{ 

StripValidBlock ( i , j , H) ; 

} 

printf ( "Stripped: %ld\r" , ValidBlockStack . BlockCurSize) ; 

} 

} 

ChipArrayClass :: St ripRandomBlocks (long Num, int H) 
{ 

long i; 



int X,Y; 



for (i=0; i<Num; i++) 

{ 

if (PickRandomValidBlock(X,Y,l,H) ) 
StripValidBlock(X,Y,H) ; 

} 

return (TRUE) ; 

1 



ChipArrayClass : : StripAreaToBlock(BlockCla3S *TempBlock, int X, int Y, int Width, int Height) 
{ 

if (X+Width>Xdim | | Y+Height>Ydim | } X<0 | j Y«:0) 
return {FALSE) ; 

// strip an area of the chip into a block 
TempBlock->Allocate (Width*Height ) ; 
TempBlock->WSize = Width; 
TempBlocfc- >HSize = Height; 

int counter = 0; 

int i , j ; 

for (i=0; i<Width; i++) 
for (j=0; j<Height; 
{ 

TempBlock->DataStack [counter] = DataGrid [X+i] [Y+j ] ; 
DataGrid[X+i] [Y+j] = NULL; // removed 
TetnpBlock- >DataStack [counter] ->SetRelative (i, j) ; 
counter++; 

} 

} 

ChipArrayClass ; .-CheckBlockFitToArea (BlockClass *TenipBlock, int X, int Y) 
{ 

int valid = TRUE; 
int i ; 
int tx, ty; 

if (TempBlock==NULL) 
return {FALSE) ; 

for (i=0; i<TempBlock->DataStackSize && valid; i++) 

{ 

if (TempBlock->DataStack [i] !=NULL) 

{ 

tx = X+TempBlock->DataStack [i] ->relx; 
ty = Y+TempBlock->DataStack [i] ->rely; 
if (tx<Xdira && ty<Ydini && tx>-l && ty>-l) 

if (DataGrid [tx] [ty] !=NULL) 
valid = FALSE; 

else 

valid - TRUE; 

else 

valid = FALSE; 

} 

} 

return (valid) ; 

} 

ChipArrayClass; :PutBlockInArea(BlockClass *TerapBlock, int X, int Y) 
{ 

int i ; 
int tx, ty; 

if { !CheckBlockFitToArea{TempBlock, X, Y) ) 
return (FALSE) ; 

for (i=0; i<TerapBlock->DataStackSize; i++) 

{ 

if (TeTnpBlock->DataStackti] !=NULL} 
{ 

tx = X+TempBlock->DataStack [i] ->relx; 



} 



ty = Y+TerapBlock->DataStack [i] ->rely ; 

DataGridftxJ ttyj = TempBlock->DataStack [ij ; 
TempBlock->DataStack [i] =NULL; 



TempBlock->DeAllocate() ; // toast! 
return (TRUE) ; 



} 

ChipArrayClass : : ChipArrayClass ( ) 

{ 

DataGrid = NULL; 
Xdira = Ydim = SynthSteps = 0; 
Radius = 9; 

MaxAl lowed = 4; 
GlobalHeight = 2; 

ScanRadius = 1 ; 
LeakageHalf Life = 1; 
weight flag = 0; 

1 

ChipArrayClass : : -ChipArrayClass ( ) 
{ 

_ DeAllocateO ; 

^ ChipArrayClass :: Allocate {int X, int Y) 

^ DataGrid = new LocalDataClass ** [X] ; 

L, if (DataGrid==NULL) 

return (FALSE) ; 
Q int i , j ; 

Z'. for (i=0; i<X; i + + ) 

^ { 

DataGrid ti] = new LocalDataClass * [Y] ; 
=j if (DataGrid [i] ==NULL) 

^ return (FALSE) ; 

Ai for (3=0; ]<Y; 

^ DataGrid [i] [j] = new LocalDataClass; // featherweight objects 

if (DataGrid [i] [j]==NULL) 
^ return (FALSE) ; 

^ } 

Xdim = X; 
Ydim = Y; 
return (TRUE) ; 



ChipArrayClass: :DeAllocate ( ) 
{ 

if (DataGrid==in;LL) 
return (TRUE) ; 
int i , j ; 

for {i=0; i<Xdim; i++) 

{ 

for (j=0; j<Ydira && DataGrid [i] f =NULL; 
{ 

if (DataGrid [i] [j] !=NULL) 

delete DataGrid [i] [j]; 

} 

if (DataGrid fi] !=NULL) 

delete[] DataGrid[i]; 

} 

delete I] DataGrid; 
DataGrid = NULL; 
Xdim = Ydim = 0; 
SynthSteps = 0; 



ChipArrayClass r :ReadCdl {char *FileName, int realflag) 
{ 

PILE *ifp; 
int maxX, tnaxY, X, Y; 
char datastring [MXLINE] ; 
int flag=TRUE; 

if (realflag) 

flag = ReadCdl (FileName, FALSE); 
if (iflag) 

return (FALSE) ; 

ifp = f open (FileName, "rt"); 
if (NULL==ifp) 
{ 

pr int f ("Unable to open: %s\n", FileName); 
exit (1) ; 

} 

fgetsidatastring, MXLINE, ifp) ; 
maxX = 0; 
maxY = 0; 

while (!feof(ifp) ! f error (ifp) ) 
{ 

fgets (datastring, MXLINE, ifp); 
if (fecf(ifp) II terror (ifp)) 
break; 

sscanf (datastring, "%d %d" , ScX, &Y) ; 
if (X>maxX) 

niaxX=X; 
if (Y>maxY) 

maxY=Y; 
if (realflag) 

DataGrid[X] [Y] ->cdldata . LineScan (datastring) ; 
if (X==0) 

printf ("%d\r" , Y) ; 

} 

fclose(ifp) ; 
if (! realflag) 

{ 

maxX++; 
maxY++; 

flag = Allocate (maxX, maxY) ; 
return (flag) ; 

} 

return (TRUE) ; 

} 

ChipArrayClass : :DumpCdl (char * FileName) 
{ 

FILE *fp; 
int i , j ; 

fp = fopen( FileName, "wt"); 
fprintf (fp, 

"X\tY\tPROBE\tDESTYPE\tFEATURE\tQUALIFIER\tEXPOS\tTBASE\tENDPOS\tPOSITION\tPBASE\tFINISHPOS\tFIXED\t 
VARIABLE\tUNIT\tBLOCK\tATOM\tREPEAT\tSEQNO\tLAYOUT\tACCESSION\tLOCUS\n") ; 
for (j=0; j<Ydim; j++) 

{ 

for (i=0; i<Xditn; i++) 

{ 

if (DataGrid[i] {j] t=NULL) 

DataGrid[i] [ j ] ->cdldata.DumpLine (fp, i,j); 

else 
{ 

printf ("Null value in grid %d %d\n", i,j); 

} 

} 

printf ("Out Cdl: %d\r" , j); 

} 

f close ( fp) ; 



ChipArrayClass : :GenerateMut File (char *FileName) 
{ 

FILE *fp; 
int i , j ; 

fp = fopen{FileName, "wt"); 
for (3=0; j<Ydim; 

{ 

for (i=0; i<Xdim,- i++> 

{ 

if (DataGrid[i] [j] !=NULL} 

{ 

DataGrid [i] ( j J ->cdldata .DumpMut (fp) ; // single descriptive character 

} 

else 
{ 

fprintf(fp, "-"); 

printf("Null value in grid %d %d\n" , i,j); 

} 

} 

fprintf(fp, "Nn"); 
printfC'MUT: %d\r", j); 

) 

f close ( fp) ; 



ChipArrayClass; :GenerateDiff Pile (char *FileName) 

( 

FILE *fp; 
int i , j ; 

int tn, te, ts, tw; 
double n,e,s,w; 
double count; 

fp = fopen{FileName, "wt") ; 

for (j=0; j<Ydim; 

{ 

for {i=0; i<Xdim; i++) 
{ 

if (Valid(i,j}) 

{ 

fprintf(fp, "X:%d\tY:%d\t", i,j); 
tn=ts=tw=te =0; 
if (Valid(i, j-1) ) 

tn - DataGrid[il [j-1] ->retdata.Diff {DataGrid[i] [ j ] - >retdata) ; 
if (Validd, j+1) } 

ts = DataGridfi] fj+lj ->retdata .Diff (DataGrid [iHjJ ->retdata) ; 
if {Valid(i-1, j) ) 

tw = DataGrid[i-l] [j] ->retdata.Diff (DataGridti] [ j ] ->retdata) ; 
if {Valid(i+1, j) ) 

te = DataGridfi+l] fj] ->ret:data.Diff (DataGrid [ij [j J ->retdata) ; 
fprintf (fp, "N:%d\tE:%d\tS:%d\tW:%d\tT: %d\t" , tn, te, ts, tw, tn+te+ts+tw) ; 

fprintf{fp, "LAST:%d\t", DataGrid [i] [j ] ->retdata .GetLast {)) ; 
fprintf (fp, "BREADTH ; %d\t" , DataGrid[il [j] ->retdata. GetLast {) -DataGrid[il [j] - 
>retdata.GetFirst 0 } ; 

DataGrid [i} [j] ->PrintLongSeq(fp) ; 

fprintf{fp, "\t%s\n", DataGrid [i] [j ]- >cdldata . qualifier) , • 

count++; 
n+=tn; 
3+=tS; 
e+=te; 

W+=tW; 

} 

else 
{ 

printfC'Null value in grid %d %d\n'<, i,j); 

} 

} 



print f {"Diff : %d %lf %lf %lf %lf %lf\r", j, n/count, e/count , s/count , w/ count, 
(n+e+s+w) / (4*count) ) ; 
} 

//fprintf (fp, "Diff : %d %lf %lf %lf %lf %lf\n", j, n/count, e/count, s/count , w/count, 
(n+e+s+w) / {4*count ) ) ; 
f close (fp) ; 

} 

ChipArrayClass : :ReadRet (char *FileName, int realflag) 
{ 

FILE *fp; 
int i , j , k ; 

char dataline [MXLINE] ; 
long total; 

if (realflag) 

ReadRe t ( F i 1 eName , 0 ) ; 
fp = fopen(Fil eName, "rt"); 
if (fp==NULL) 
{ 

print f ( "Unable to open: %s\n", FileName) ; 
exit (1) ; 

} 

if (realflag) 
{ 

for (i=0; i<xdi[n; i++) 
for (j=0; j<Ydim; 
{ 

if (DataGridli] [j] !=NULL) 

DataGrid[i] [j ]- >retdata .Allocate (SynthSteps) ; // allocate this data 

else 

printf { "Death by lack of allocation\n" ) ; 

} 

} 

k=-l; 

j=Ydim; 
total = 0; 

while (fgets (dataline, MXLINE, f p) ) 
{ 

if (dataline [0] =='r' ) 
{ 

sscanf (dataline, "reticle: %s", &retname) ; // set up reticle template name 
retname [strlen(retname) -2] = '\0'; 

// initialize for reading next lines 

k++; 

printf ("Reticle: %d\r", k) ; 
j=Ydim; 
continue; 

} 

if (strlen (dataline) <10 1| dataline to] ==';' ) 
continue; 

if (dataline [0]=='0' || dataline [0] ==' 1 ' ) 



{ 



3--; 

for (i=0; icXdim 3>-l; i++) 

if (dataline [i] ==' 1' ) 
{ 

if (realflag) 

{ 

DataGridfi] [j] ->retdata. SetBit (1, k) ; 

} 

total ++; 



NumOnes = total; 
SynthSteps = k+1; 

f close (fp) ; 



ChipArrayClass : : DumpRet { char *FileNarae) 
{ 

FILE *fp; 

int 

char dateline [POCLINE] ; 

Ep = fopen (FileName, "wt") ; 

fprintf (fp, This file has been annealed to minimize edges\n" 

for (k=0; k<SynthSteps; k++) 
{ 

fprintf (fp, "\n\n;base: X") ; 

fprintf (fp, "\nreticle: %s%02d", retname, (k+l)); 
fprintf (fp, "\nR I 1 1 0 0 %d %d %d",Xdim, Ydim, 1) ; 
for (j=Ydira-l; j>-l; j--) 
{ 

fprintf (fp, "\n"); 

for (i=0; i<Xdim; i++) 

{ 

if (DataGridfi] [j]==NULL) 
{ 

prlntfC'Data leakage: %d %d\n", i,j); 

} 

else 

if (DataGrid[i] [j] ->retdata.GetBit (k) ) 
dataline[i] = '1'; 

else 

datalinefij = '0'; 

} 

dataline [Xdim] = ' \0 ' ; 
fprintf (fp, "%s", dataline); 

} 

fprintf (fp, "; \nO;\n") ; 
print f ("DUMPRET: %d\r", k) ; 

} 

f close (fp) ; 

} 

ChipArrayClass: : InterpretlnstructionLine (char *Line) 
{ 

char TempStr [MXLINE] ; 
int height; 
long searchlimit; 
long start, finish; 
int tx,ty,x,y; 
int value; 
int destype; 
int radius, max; 
dovible dval ; 

// read an instruction and do the appropriate thing 
if (Line[0]==' ; ') 

return (TRUE) ; // comment 
3scanf{Line, "%s", TempStr); // pick off the initial piece 
if {! strcmp (TempStr, "READCDL:")) 
{ 

Sscanf(Line, "READCDL: %s", TempStr); 
ReadCdl (TempStr, TRUE) ; 
return (TRUE) ; 

} 

if (! strcmp (TempStr , "READRET:")) 
{ 

s3canf(Line, "READRET: %3", TempStr); 
ReadRet (TempStr, 1) ; 



• 



i» 



return (TRUE) ; 

f { !strcrnp(TempStr, "DUMPCDL:")) 

sscanf(Line, "DUMPCDL: %s", TempStr} ; 
Dumped 1 (TempStr) ; 
return (TRUE) ; 

f { 1st rcrap (TempStr, "DUMPRET:")) 

sscanf(Line, "DUMPRET: %s", TempStr); 
DumpRet (TempStr) ; 
return (TRUE) ; 

f ( 1st rcrap (TempStr, "DUHPMUT:")) 

sscanf(Line, "DUMPMUT: %s", TempStr) ; 
GenerateMut File (TempStr) ; 
return (TRUE) ; 

f ( (strcmp (TempStr, "DUMPDIFF: ") ) 

sscanf(Line, "DUMPDIFF: %s", TempStr); 
GenerateDiff File (TempStr) ; 
return (TRUE) ; 

if (! St rcmp (TempStr, "STRIPBLOCKS : " ) ) 

sscanf(Line, "STRIPBLOCKS: %d" , &height) ; 
GlobalHeight = height ; 
St ripAllValidBlocJcs (height) ; 
Shuffle 0 ; 
return (TRUE) ; 

f ( ! Strcmp (TempStr, "DIAGONALRE PLACEMENT : " ) ) 

sscanf(Line, " DIAGONALRE PLACEMENT : %ld" , &searchlimit) ; 
DiagonalReplacement (searchlimit ) ; // put all bloclts back on chip 
return (TRUE) ; 

f {! strcmp (TempStr, "HORIZONTALREPLACEMENT : " ) ) 

sscanf(Line, "HORIZONTALREPLACEMENT: %ld" , &searchl imit ) ; 
DiagonalReplacement (searchlimit ) ; // put all blocks back on chip 
return (TRUE) ; 

.f ( ! strcmp (TempStr, "AGGREPLACEMENT : ") } 

sscanf(Line, "AGGREPLACEMENT: %ld" , &searchlimit) ; 
AggregateReplacement ( searchlimit) ; // put all blocks back on chip 
return (TRUE) ; 

f {! strcmp (TempStr, "SETVALIDUNITS : " ) ) 

sscanf(Line, "SETVALIDUNITS: %ld %ld %d" , &start, ^finish, &value) ; 
SetUnits (start, finish, value); 
return (TRUE) ; 

f ( ! strcmp (TempStr, "SETVALIDAREA: ") ) 

sscanf(Line, "SETVALIDAREA: %d %d %d %d %d" , &x, &y, &tx, &ty , &value) ; 
SetArea (x, y, tx, ty, value) ; 
return (TRUE) ; 

,f (! strcmp (TempStr, "SETVALIDANTIT^EA: " ) ) 

sscanf (Line, "SETVALIDANTIAREA: %d %d %d %d %d" , &x, &y , &tx, &ty, lvalue) 
SetAntiArea (x, y, tx, ty. value) ; 
return (TRUE) ; 



if ( 1st rcmp (TempStr, "SETVALIDDESTYPE : " ) ) 



{ 

sscanEILine, "SETVALIDDESTYPE : %ld %d" , &destype. &value) ; 
SetDestype (destype, value); 
return (TRUE) ; 

if { !strcmp(TempStr, "STRIPBADPROXIMITY: " ) ) 

ascanf {Line, "STRIPBADPROXIMITY: %d %d %d" , &height, tradius, &max) 
Radius = radius; 
MaxAl lowed = max; 
StripBadProximity Values (height) ; . 
return (TRUE) ; 

if ( !3trcmp(TetnpStr, "SETPROXIMITY: " } ) 

sscanf(Line, "SETPROXIMITY: %d %d" , Sradius, &raax) ; 
Radius = radius ; 
MaxAl lowed = max; 
return (TRUE) ; 

if ( ! strcTnp(TempStr, "FIXBAD:")) 

sscanf{Line, "FIXBAD: %ld", Sisearchlimit) ; 
DoubleReplacement (searchlirait) ; 
return (TRUE) ; 

if ( !3trcmp(TeinpStr, "WEIGHT:")) 

sscanf(Line, "WEIGHT: %d%lf", fcradius, &dval) ; 
ScanRadius = radius; 
LeakageHalf Lif e = dval; 
weightflag = TRUE; // use weights 
return (TRUE) ; 

if ( ! strcnip(TempStr, "NOWEIGHT: ") ) 

weightflag = FALSE; 
return (TRUE) ; 

return (FALSE) ; 

} 

ChipArrayClass : : Readlnstruct ionFile (char *FileName] 
{ 

// read the file and be happy 
FILE *fp; 

char dataline[MXLINE] ; 

fp = f open (FileWame, "rt"); 
if (fp==NULL) 

return (FALSE) ; 
while (fgets (dataline, MXLINE, fp) ) 
{ 

InterpretlnstructionLine (dataline) ; 

} 

f close (fp) ; 



TestTwo{ ) 

ChipArrayClass Test; 
Test .ReadlnstructionFile ( "test .opt ") ; 



Live (char *FileName) 



ChipArrayClass Test; 
Test .ReadlnstructionFile (FileName) ; 



main(int argc, char **argv) 
{ 

if {argc==2) 
{ 

Live <argv [1] ) ; 

} 

else 

{ 

printfC'File name required 

} 



